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ABSTRACT 

We present a study of close white dwarf and M dwarf (WD+dM) binary systems and examine the effect 
that a close companion has on the magnetic field generation in M dwarfs. We use a base sample of 1602 
white dwarf - main sequence binaries from Rebjissa-Mansergas et aTJ ( 120101) to develop a set of color cuts 
in GALEX, SDSS, UKIDSS, and 2MASS color space. Then using the SDSS DR8 spectroscopic database, 
we construct a sample of 1756 WD+dM high-quality pairs from our color cuts and previous catalogs. We 
separate the individual WD and dM from each spectrum using an iterative technique that compares the WD 
and dM components to best-fit templates. Using the absolute height above the Galactic plane as a proxy for 
age, and the Ha emission line as an indicator for magnetic activity, we investigate the age-activity relation for 
our sample for spectral types < M7. Our results show that early-type M dwarfs (< M4) in close binary systems 
are more likely to be active and have longer activity lifetimes compared to their field counterparts. However, 
at a spectral type of M5 (just past the onset of full convection in M dwarfs), the activity fraction and lifetimes 
of WD+dM binary systems becomes more comparable to that of the field M dwarfs. One of the implications 
of having a close binary companion is presumed to be increased stellar rotation through disk-disruption, tidal 
effects, or angular momentum exchange. Thus, we interpet the similarity in activity behavior between late-type 
dMs in WD+dM pairs and late-type field dMs to be due to a decrease in sensitivity in close binary companions 
(or stellar rotation), which has implications for the nature of magnetic activity in fully-convective stars. Using 
the WD components of the pairs, we find WD cooling ages to use as an additional constraint on the age- 
activity relation for our sample. We find that, on average, active early-type dMs tend to be younger and that 
active late-type dMs span a much broader age regime making them indistinguishable from the inactive late-type 
population. We also show that magnetic strength, as measured by Ha, is comparable between paired and field 
M dwarfs until a spectral type of M6/M7 where M dwarf activity for stars with close companions becomes 
much stronger. In addition, we present 37 very-close candidate pairs with fast-moving orbits that display radial 
velocity changes over hour time-scales. 

Subject headings: Stars: low-mass — Stars: white dwarfs — Stars: activity — Stars: rotation — binaries: 
spectroscopic — binaries: close 



1. INTRODUCTION 

Aside from being cooler, dimmer, and much less mas- 
sive than the stars like the Sun, M dwarfs (dM) are the ma- 
jor stellar constituents of the Galaxy. The ubiquitous nature 
of dMs, in addition to their long lifetimes, make them ex- 
celle nt tracers of Galactic structure, kin ematics, and evolu- 
tion (iGiiisetiD USqI L^HeFaD |20Ql). In addition, dMs 
have strong magnetic fields that heat their upper atmospheres, 
producing excess emission in the chromosphere and corona 
(called "magnetic activity"). This magnetic activity plays 
a vital role in the space weather environments of dMs and 
likely affects the true habitability of attending planets. How 
these small stars generate such strong fields (^several kG) 
and the subsequent magnetic activity is still not well under- 
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stood. Therefore, detailed observations of the magnetic prop- 
erties of dMs are vital for constraining theoretical dynamo 
models. Adding to their utility, the dM sequence (M0-M9) 
also probes an important transition in the structure of stellar 
interiors, where they change from having both radiative and 
convective zones, to being fully convective (around a spectral 
type of M4: e.g.. iDantona & Mazzitellil 1 1 985t iDorman et alj 
Il989t IChabrier & Baraffelll997l) . This transition has impor- 
tant ramifications for how magnetic fields are g enerated and 
maintained across the dM stellar sequence (e.g . , iMorin et"aT] 
l2008tlReiners & Basrill2009l: iMorin et al.ll2010b . 

In the Sun, magnetic field generation is thought to be pri- 
marily driven by the shear between the stratified solid body 
rotation of the radiat ive zone and the differential rotation o f 
the convective zone (lParkerlll955L[l99l lO ssen drijverll200l . 
This transition region, called the tachocline, acts to strengthen 
and stretch poloidal magnetic fields, establishing a large-scale 
toroidal magnetic field that can become unstable and rise 
through the convective zon e to become an active region on 
the surface of the star (e.g.. lBrowni ng 2008). Because of the 
importance of the tachocline shearing in this model, stellar ro- 
tation plays a major role in the strength of the magnetic field 
and resulting activity. Magnetic field generation in early-type 
dMs (< M4) is also thought to be due to a similar mechanism 
as early-type dMs still have tachoclines. However, dMs later 
than ~M4 (late-type) are modeled to be fully convective (e.g., 
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IChabrier & Bara ffe 1997) and lack tachoclines. Despite these 
structural differences, many studies have shown that late-type 
dMs have strong fields and are magnetically ac tive (the frac- 
tion of stars that are active peaks around M7; [Hawlev et all 
[l996t iGizITet all 120001: IWest et al.N200l 1200a 120081 here- 
after W08~ TDonati et al.ll2008t IReiners & B asrj 120071 120091 
l2Ol0l) 

In most stars, the strength (and presence) of magnetic 
activity appears to correlate with age a nd rotation rate (as 
sugges ted above). Several decades ago, IWilson & WoollevI 
d 19701) discovered a relation linking magnetic activity 
(emission in Ca II H&K) and the age of a star (estimated 
from galactic orbital parameters) using main sequence stars 
(mostly G dwarfs). They found that stars that had undergone 
many dynamical interactions (i.e., are older) had lower levels 
of activity than those stars that had more circular orbits 
(i.e., are younger). A short time later, ISkumanichl (119721) 
discovered an age-rotation-activity relation (again using 
mostly G dwarfs) showing that stars tend to spin down as they 
age as well as decrease in activity. The Skumanich relation 
occurs because angular momentum is lost from magnetized 
stellar winds; as a result, magnetic fields (and activity) are 
reduced in magnitu de. Subsequent studies have confirmed the 
ISkumanichl d 19721) resul ts and establis h ed relations between 
an d activity (e.g., IWielenl 119771: iGiampapa & Liebertl 
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2008 ) as well as a ge and rotation (e.g., iBarrvl 1988 



Soderblomet al. 1991; Mohantv & Basil 2003; Barnes 2007; 
Kiraga & Stepienl 120071: IMamaiek & Hillenbrand! I2OO8I: 
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and early-type M dwarfs. 

The exact mechanism behind the magnetic field gener- 
ation in late-type dMs is still unclear. A recent study 
( IChabrier & Kiik er 2006) suggested that magnetic activity in 
late-type dMs is due to an a 2 mechanism, where flux tubes are 
twiste d in a stratified ro tating fluid acted upon by the Coriolis 
force. iBrownind (12008b used simulations to show that strong 
magnetic fields are generated in fully convective envelopes 
when large slowly overturning flows are acted upon by Cori- 
olis forces. It appears that no matter the mechanism, stellar 
rotation still plays an important role; rotation has a strong in- 
fluence on the magnitude of the Coriolis force acting upon the 
overturning flows. 

There is some observational evidence that the age-activity 
and rotation -activity relations ex t ends into the la t e-type 
dM regime dDelfosse et all 119981: iMohantv & Basril 120031: 
IReiners & Basrill2007ll2008t W08). However, there are some 
clear difference between early and late-type dMs. W08 cali- 
brated a "Galactic stratigraphy" technique (stars farther from 
the Galactic mid-plane are statistically older) to show that ac- 
tivity decreases as a function of age in all dM spectral types 
and that the activity lifetimes of M dwarfs change across the 
transition to fully con vective interiors (f r om ~ l-2 Gyr to ~7- 
8 Gyr). In addition, IReiners & B asri (2008) demonstrated 
that the vast majority of late-type dMs are rotating quickly 
(as compared to a distribution of rotation rates for early-type 
dwarfs). These results suggest that if rotation is a strong 
driver of magnetic activity, then the spindown times of late- 
type dMs are much longer than those of more massive stars. 
In a few cases, there have been reports of fast rotating dMs 
that do not show measurab le magnetic activity (IWest & Basril 
[20091 IReiners et al.ll2012b . While some of these result s may 
be due to spurious measurements (IReiners et al.ll2012l) . they 



leave open the question of exactly what role rotation plays in 
the generation of magnetic field in fully convective stars. 

M dwarfs with close binary companions may shed some 
light on the role that rotation plays in the generation, 
strength and lifetime of magnetic activity in dMs. It is 
thought that the rotational evolution of isolated stars is 
intric ately tied to the lifetime of the circumstellar disk 
(e^ lBarnes & Sofidfl996t IBouvier et alii 19971: iTinker et al.l 
2002). Longer lived disks allow for the star to lose 
more angular mome n tum through m agnetic disk-locking 
(e.g., iKoenigll 119911 : I Shu et all 119941) or through accre- 
tion powered stellar winds (e.g., [Korvcanskv & Papaloizoul 



119951: iPapaloizou & Terquemlll997l: iMatt et al.l 12012b : stars 
with long-lived disks tend t o be slower rotato rs than 
stars with short-lived disks dMeibom et al.l 120071) . Bi- 
nary compani ons at < 100 AU shorten the disk dissipa- 
tion timescale (|Lin et al.l | 1993t lArtvmowicz & LubowHl994t 
lArmitage & Clarke! 119961) . causing angular momentum loss 
mechanisms such as accretion powered stellar winds and 
magnetic disk-locking will be halted, allowing for stars in 
close binary pairs to remain rotating faster for longer. Close 
binary pairs with separations <10 AU will also undergo the 
effects of tidal forces, and over time wil l have their spin rate 
synch ronized to the orbital motion (e.g.. lMeibom & Mathieul 
120051) . increasing the stellar rotation. Therefore, stars in bi- 
nary systems with separations <100 AU should rotate faster 
due to tidal effects and/or disk truncation mechanisms. If 
rotation indeed plays an important role in the production of 
magnetic acitivity of low-mass stars, then dMs in close binary 
pairs should be more active than their single star counterparts. 

Some previous studies have found evidence for in- 
creased magnetic ac tivity in close, eclipsing dM binaries 
dMorales et alj|2010l and references therein). However, the 
paucity of low-mass eclipsing systems precludes any statis- 
tically significant study that focuses on how a close com- 
panion affects the observable properties of dMs. Due to the 
intrinsically faint nature of dMs, large samples of both sin- 
gle or binary stars were historically untenable. However, 
with the advent of larg e, all-sky surveys (e.g ., Sloan Digi- 
tal Sky Survey (SDSS, lAbazaiian etaITl2009l) Two Micron 
All Sky Survey (2MASS7 ISkrutskie et al.ll2006b . UKIRT In- 
frared Deep Sky Survey (UKIDSS^), NASA's Galaxy Evolu - 
tion Explorer (GALEX, iMartin & The GALEX Teamll2005l) . 
etc.), large spectroscopic and photometric catalogs of low- 
mass stars have emerged. SDSS alone has enabled cata- 
logs that contain over 70,000 visually confirmed dM spectra 
dWestetalJl2011l hereafter Wl l) and more than 30 million 
5-band images (LBochanski et al. 2010), providing statistically 
robust samples of field dMs for Galactic and stellar studies. 
While large samples of wide binary com panions have been 
discovered in SDSS (iDhital et al.l 120101) , identifying large 
samples of close binary companions without high-resolution 
spectroscopy or adaptive optics imaging is challenging. Sev- 
eral SDSS studies have identified samples of close dM pairs, 
but have not examined t heir effects on the measured magneti c 
activity properties (e.g jBecker et al.1201 UlClark et al.l2012l) . 

One technique for identifying close companions to dMs 
is to focus on systems containing an dM and a white dwarf 



8 The UKIDSS project is defined in Lawrence e tlul <2007l). UKIDSS uses 
the UKIRT Wide Field C amera (WFC AM Casali et al. 2007) and a photo- 
metric system described in Hewett et al. 1 20(3). The p ipeline processing an d 
science archive are described in Irwin et al (in prep) and[Hambly et all 12003) . 
We have used data from the 2nd data release, which is described in detail in 
Wan-en et al. (2007). 
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(WD), which can be easily identified in low-resolution spec- 
troscopy. Previous SDSS spectroscopic studies used close and 
unresolved WD+dM binary pairs to investig ate the magnetic 
activity of dMs with close companions ( e. g.. IRavmond et ail 
12001 ISilvestri et al1l2006l hereafter S06; iHeller et al.1120091 
S06 compared the activity fraction of 747 close WD+dM 
binary pairs as function of spectral type to that of isolated 
dMs (W04). In comparison, S06 found that the activity in 
early-type dMs is higher in cl ose binary pairs tha n in isolated 
dMs; similar to the result of IHeller et al.l (120091) who found 
23.2% of early-type unresolved white dwarf-main sequence 
(WDMS) to be active comp ared to ^10% of active early-type 
dMs from lWest et al.l ( 120041) . From the WD cooling ages, S06 
found that all of the active, early-type dM were found around 
young WDs. There was some evidence that late-type dMs 
(in WD+dM pairs) also had a higher activity fractions, but 
there were not enough p airs to statisticall y confi rm this result. 
In addition, neither the IRavmond et al.l (12003b nor the S06 
samples were large enough to put the WD+dM pairs in their 
proper Galactic context (how the magnetic activity varies with 
location in the Galaxy). Additional SDSS data releases have 
dramatically increased the sample si zes of WD+dM pairs 
dRebassa-Mansergas et al.ll2010l 120121) , which now allow for 
statistical investigations of the magnetic activity properties of 
low-mass stars in close binaries. 

In this paper, we use a large sample of WD+dMs to build 
on the S06 analysis and investigate the magnetic activity prop- 
erties of close WD+dM pairs as a function of spectral type, 
Galactic height (a proxy for age), WD cooling ages, and bi- 
nary separation. In Section 2 we describe our sample selection 
and general properties of our sample. In Section 3 we describe 
the procedures we developed to separate the binary compo- 
nents and calculate WD+dM parameters. We also present a 
detailed analysis of the WD spectra using cooling tracks to 
estimate age, mass, and temperature. In addition, Section 3 
details a radial velocity (RV) analysis using individual SDSS 
exposures to uncover a population of very close binary pairs. 
In Section 4, we present our results on activity fractions, life- 
times, and strength as a function of binary, stellar and Galactic 
parameters. Section 5 discusses the results from our analysis 
and the potential implications of having close companions on 
the magnetic activity of dMs. Finally, we summarize our find- 
ings in Section 6. 

2. DATA 

We identified 1756, high-quality, close WD+dM spectro- 
scopic binaries from the SDSS Data Release 8 (DR8). Our 
goal was not to produce a complete sample of WD+dM 
binary pairs (something that SDSS spectral selection pre- 
cludes), but rather produce a sample of high fidelity and 
high signal-to-noise (S/N) spectroscopic WD+dM binary 
pairs, expanding upon previous close W D+dM catalogs: the 
IRebassa-Mansergas et all (120101 120121) sample of 1602 and 
2248 WDMS pairs, respectively; and the S06 catalog of 747 
close WD+dM pairs. In particular, we were interested in in- 
creasing the number of late-type dMs (M7-M9) found by pre- 
vious studies in an attempt to better understand how activity in 
close pairs affects fully convective M dwarfs at the bottom of 
the main sequence. We employed the use of multi-wavelength 
color-cuts (ultraviolet, optical and infrared) to produce a sam- 
ple of 1756 close WD+d M spectroscopic binarie s selected 
from the SDSS DR8 (DR8. lAbazaiian et al.l l2009l) using the 
IRebassa-Mansergas et alj ( 120101) WDMS sample as a guide- 
line for developing our own color cuts (described below). 



The advantage of using a photometric selection criteria that 
spanned from ultraviolet (UV) to infrared (IR) wavelengths 
was an increased sensitivity to late-type pairs at infrared 
wavelengths, whose optical fluxes are dominated by the WD 
components. Inversely, this multi-wavelength technique was 
also sensitive to pairs with very faint WD components whose 
optical flux is dominated by the dM. In the end, our sample of 
17 56 close WD+dM pairs does not contain as many pairs as 
the IRebassa-Mansergas et al.l (120121) sample due to our strin- 
gent quality cuts that ensure high S/N for both the WD and 
dM component; which were necessary for our magnetic ac- 
tivity analysis. 

We used photometry from four different public catalogs, 
spanning the ultraviolet, optical, and infrared wavelengths, 
in addition to the optical spectra from SDSS DR8 cata- 
log. The/ M v (1350-1750 A) and nuv (1750-2800 A) bands 
from GALEX; the uqriz bands (3000- 1 1000 A) from SDSS 
dFukugita et al.lll996tlGunn et al.lll998l) : and the J (1.25/im), 
H (1.65/im), and K (2.15/xm) band from both 2MASS and 
UKIDSS were used. I n addition to sl i ghtly different filters, 
2MA SS (Vega-based: ISkmtskie et all l2006h and UKIDSS 
(AB; iHewett et alj 1200617 use different magnitude systems. 
While 2MASS is an all-sky surey, UKIDSS probes a much 
smaller area but goes three magnitudes deeper, creating an 
effective volume 12 times larger than 2MASS. Therefore, 
using both 2MASS and UKIDSS maximizes the number of 
WD+dMs that are detected in the infrared. 

2.1. Sample Selection 

We selected all SDSS "point sources" with spectra that sat- 
isfied the following quality cuts: ugriz < 24 and redshift < 
0.003. Our initial leniency helped facilitate a large number 
of candidate pairs, while leaving room to filter out low qual- 
ity objects in later stages of the sample selection. After re- 
moving duplicate objects, our initial sample was composed of 
520,701 SDSS stars. Multi-wavelength counterparts to these 
objects were found in GALEX, 2MASS, and UKIDSS us- 
ing three arcsecond matches to their positions. There were 
227,728 matches in GALEX, 300,665 matches in 2MASS, 
and 1 3 1 ,4 1 6 matches in UKIDS S . 

We reduced the number of false positive candidate 
WD+dM binary systems by enforcing additional color-cuts 
in fuv, nuv, u, g, r, i, z, J, H, and K color space. We used 
the color cuts provided by Reba ssa-Mansergas et al.l (120101) 
as an initial guide for our se lection. We visually inspected the 
1602 WDMS spectra from IRebassa-Mansergas et al.l (120101) 
and concluded that 1399 were WD+dM pairs. We then ex- 
plored the color space defined by all of the above filters and 
traced out regions with minimimal contamination (we de- 
fined contamination to mean the percentage of non WD+dM 
pairs included in a color cut). Our final color cuts were 
simila r to those described in the Rebassa-Mansergas et al. 
(120101) analysis, but were more useful in probing later spec- 
tral types (M7-M9) as well as spectra dominated by WD or 
dM flux. Our color cuts were applied to the latest SDSS 
DR8 catal og, where we found 364 new p airs that were not 
present in IRebassa-Mansergas et al.l (120101) . For the follow- 
ing discus sion we will compare the effectiv eness of our color 
cuts to the lRebassa-Mansergas et aD (120101) sample (on which 
we base our sample selection). Since we began our anal- 
ysis, a larger sample of WD+d M pairs has been published 
( IRebassa-Mansergas et al.ll2012l). We will com pare our sam- 
ple to that of IRebassa-Mansergas et all ( 120121) in a later sec- 
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Figure 1. u - z versus r — z color-color diagram for stars in 
SPSS. The colored filled circ les represent known WD+dM pairs from 
Rebassa-Mansergas et al. 1 2010). The symbols are colored according to the 
spectral type of the M dwarf. The color cut is represented by the black- 
dashed line (all objects below that line are selected). The grey contours are 
all of the 520,701 stars in the initial SDSS sample; the contour levels are 10, 
100, 1000, 5000, and 10000 objects. Since all of our objects have SDSS col- 
ors, these colors provide the most comprehensive cut. We choose r - z as its 
sensitive to dM spectral ty pe while u — g isola tes the WD+dM population 
from the main stellar locus (Smolcic et al. 2004). 



tion. The color-cuts are shown in Figures HH In each figure, 
the contoured data represent all SDSS point sources, while 
the colored filled circles represe nt confirmed WD+dM from 
dRebassa-Mansergas et al.ll2010l) . The circles are colored ac- 
cording to their spectral type; this was done in an effort to 
identify trends in the spectral types in order to best select re- 
gions hosting the later M7-M9 spectral types that have not 
been well sampled in previous studies. 

When using J, H, an d K colors, we ne eded to distin- 
guish between 2MAS S dHewett et alj 120061) and UKIDSS 
(ISkrutskie et al.ll2006h photometry separately due to the dif- 
ferences in the photometric precision, the zeropoints, and the 
filter responses. Treating them separately did not identify any 
significant trends that could be due to using slightly different 
photometric systems. However, we will present our 2MASS 
and UKIDSS color cut analysis separately below. 

2.2. Color Cuts 

WD+dM systems are well separated from main sequence 
stars in a r — z vs. u — z color-color diagram (ISmolcic et al.l 
120041) . This is largely due to the fact that the u — z color 
spans the wavelengths where the WD and dM components 
peak. The r — z color serves to segregate the M dwarfs as a 
function of spectral type. Figure Q] shows the r — z vs. u — z 
color-color diagram for all SDSS stars. The main-sequence 
stars are shown as grey contours, while the WD+dM sys- 
tems s pectroscopically identified in Rebassa-Mansergas et al 
(l2OT0l) are shown as circles, with the colors corresponding to 
spectral type of the dM. The WD+dM population is clearly 
segregated from the main-sequence stars in u — z. We did a 
by-eye fit to separate the WD+dM pairs from the main stellar 
locus, shown as dashed lines in Figure[T]and described by: 

( 0.767 (r-z) + 0.786; 

for (r-z) < 0.416, 
4.573 (r-z)- 0.798; 

for 0.416 <{r-z)< 1.117, 
0.888 {r-z) + 3.318; 

for (r-z)> 1.117. 



(it - z) < < 



(1) 



As can be seen in Figure Q] this selects the vast majority of 
WD+dM pairs. However, systems that are dominated by ei- 
ther the WD or dM component might not be separated in the 
optical wavelengths. Using the UV (where the WD domi- 
nates) and the near-IR (where the dM dominates) separates 
those WD+dM systems more cleanly. Therefore, we used the 
UV data from GALEX and the near-IR 2MASS and UKIDSS 
to identify additional pairs. 

Figure [2] shows the fuv — nuv vs. J — H plot for our 
sample. The fuv — nuv isolated objects with significant UV 
flux, while J — H color separated the M dwarfs from the G 
and K dwarfs. Whil e this cut is similar to th e fuv — nuv 
vs. r — z cut used in Rebassa-Manserg as et all ( 120101) . using 
the J — H color, instead of r — z color, produced a cleaner 
separation. The locus of the WD+dM pairs was described by 



(J — -ff)2MASS > 



' 7.378 (fuv - nuv) - 1.754; 

for (fuv — nuv) < 0.243 
0.104 (fuv - nuv) + 0.013; 

for (fuv — nuv) > 0.243, j 



>(2) 



for 2MASS colors and by 

(J - i?)uKiDss > 0.068 (fuv - nuv) + 0.103 (3) 

for UKIDSS colors. We then used the nuv — r vs. nuv 
— H diagram to further separate WD+dM pairs from the 
stellar locus. As Figure [3] shows, WD+dM pairs are above 
the main-sequence stars in nuv — H while nuv — r spreads 
the pairs as a function of the dM spectral type. This cut is 
only slightly different from the nuv — i vs. nuv — H used in 
iRebassa-Mansergas et al.l (120101) and is a better tracer of dM 
spectral type. The derived WD+dM locus is described by 



#2MASS > 



for 2MASS colors and 



1.391 (nuv -r) + 0.960; 

for (nuv -r) < 2.226, 
1.147 (raw -r) + 1.513; 

for (nuv -r) > 2.226, 



(4) 



-ffuKross > 1-339 x (nuv 



0.888. (5) 



for UKIDSS colors. As the fourth, and final, color cut, we 
used nuv — K vs. z — J, which is only slightly different from 
the third one we used. However, z — J is more sensitive to sys- 
tems with later-type dM secondaries, as it includes the wave- 
length regime where late-dMs peak in flux. Figure 2] shows 
the WD+dM locus well separated above the main stellar se- 
quence. The derived WD+dM locus is described by 



z — J2MASS > 



0.198 (nuv - 
for (nuv 
0.032 (nuv - 
for (nuv 



^2MASs) + 0.220; 

- K 2 mass) < 3.813, 
#2MAss) + 0.853; 

- K 2 mass) > 3.813, 



>(6) 



for 2MASS colors and by 



Z — ./UKIDSS > 



0.175 (rati; - 
for (nuv 
0.048 (nuv - 
for (nuv 



^ukidss) + 0.388; 

- ^ukidss) < 3.457, 
i^UKiDss) + 0.827; 

- ^ukidss) > 3.457, J 



>(7) 



for UKIDSS colors. 

Using the color cuts above, we selected 53,279 unique ob- 
jects. The color cuts were able to recover 1379 of the 1399 
WD+dM pairs from Rebass a-M ansergas et al. (20lti|); the 
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Figure 2. fuv - nuv vs. J - H color-color diagram. As in Fig. 1, the filled circles are colored by M dwarf spectral type for known WD+dM pairs. On the 
left panel, the contour represent the 2MASS data and there is only one contour levels representing stellar densities >5. The right panel uses UKIDSS data and 
only uses contour levels with stellar densities >3. The distinction was made because of the photometric depth and precision disparity between the two samples. 
This cut was chosen so that the fuv - nuv colors would isolate new dM flux dominated pairs, where the WD may not component may not be obvious in optical 
photometry but will show up in the UV regime. Similarly, J - H selects WD flux dominated pairs where the dM will reveal itself in the infrared. J - H is also 
sensitive to late-dM spectral types. 
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Figure 3. nuv - r vs. nuv - H color-color diagram, filled circles are colored by spectral type as before as well as distinguishing between 2MASS and UKIDSS 
separately. The contours levels on the right panel represent stellar densities of 5, 25, 100 , 250, and 500 in the 2MASS data and 5, 20, 50, and 100 in the UKIDSS 
data on the left panel. This cut was chosen for its ability to separate the WD+dM bridge I Smolcic et al. 2004) from the main stellar locus. Colors nuv - r provides 
excellent spread in spectral type, while nuv - H helps select dMs with later spectral types. 



pairs that were missed by our color cuts had inconsistent col- 
ors with the WD+dM locus (seen in Figure [TJ. Each object 
was visually inspected and given one of three simple clas- 
sifications: not a WD+dM, a high probability candidate, or 
a low probability candidate. The first two classifications are 
self explanatory and the third was necessary for objects that 
were difficult to differentiate from single WDs or dMs with 
low S/N or other astronomical objects. After classifying each 
object by-eye, there were 1,891 high-probability WD+dM 
pairs and 2,220 low-probability pairs. We performed a de- 
tailed spectral analysis (described in Section 3) on all of the 
candidate pairs and were able to confirm that of the 1,891 
high-probability pairs 1,460 were real and of the 2,220 low- 
probability pairs 199 were real, resulting in an initial sample 
size of 1659 WD+dM pairs. The r — z vs. u — z cut se- 
lected 1605 of the WD+dM while the IR and UV color cuts 
were able to select the remaining 54 pairs that were not in- 



cluded in the r — z vs. u — z cut (Equation [T]). Of those 54 
remaining pairs, only one pair was found exclusively using 
the fuv — nuv vs. J — H cuts (Equations[2]and[3]l, 18 pairs 
were found exclusively from the nuv — r vs. nuv — H cuts 
(Equations @] and [5]), and three pairs were found exclusively 
from the nuv - K v s. z -J cuts (Equations [6] and [7]). 

In creating the IW1 ll M dwarf sample, over 70,000 dM 
were visually confirmed from over 116,000 candidate spec- 
tra from SDSS DR7. During the visual inspection process, 
candidate WD+dM objects were flagged for later analysis, 
from these we added 194 possible WD+dM that were not 
in our sample. In addition, we mat ched our sample to those 
of iRebassa-Mansergas et all (120101) and S06 and found 522 
and 144 objects, respectively, that we had missed. This re- 
sulted in a total of 910 candidate WD+dM objects that we 
had missed in our color cut analysis. After visually inspect- 
ing each spectrum, we added 141 objects to our sample. The 
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Figure 4. nuv - K vs. z — J color-color diagram. Filled circles are colored by M dwarf spectral type. Left and right panels distinguish between 2MASS and 
UKIDSS data respectively. The contours levels on the right panel represent stellar densities of 5, 25, 100, 250, and 500 in the 2MASS data and 5, 20, 50, and 100 
in the UKIDSS data on the left panel. These cuts were chosen specifically to identify late spectral types with z — J, which covers the wavelength regime where 
there is a peak in M7-M9 flux. 

other 769 objects were eliminated as they were non-WD+dM 
pairs, WD or dM stars with low S/N in the red or blue portion 
of the spectrum (making separation of a companion difficult), 
or WD+dM pairs with S/N too low for our magnetic activ- 
ity analysis. Finally, we did a final inspection of our sample 
and found 44 objects that had been misclassified or where the 
S/N was less than 3 near Ha absorption line and therefore did 
not meet our quality criteria. This resulted in a final sample of 
1,756 confirmed close WD+dM pairs. As we have mentioned 
earlier, this sample was not devised to be complete but rather a 
sample of high fidelity and high signal-to-noise spectroscopic 
WD+dM pairs. The complete catalog will be available on- 

line|I3 

When we compared our fin a l sam ple to the S06 and 
iRebassa-M ansergas et alj (120101 120121) catalogs, we found 
that we were missing 310, 409, and 606 pairs, respectively. 
However, we added 847, 556, and 121 new pairs, respectively, 
to those three samples. The pairs missing from our sample 
were visually inspected; and the majority of them were not 
WD+dM pairs but rather WD+K pairs, had low S/N, or the 
components could not be properly separated with our spec- 
tral analysis procedure (described in Section 3.1). Thus, they 
were not useful for our magnetic activity study and not in- 
cluded in the final sample. 




3. ANALYSIS 

3.1. WD+dM separation procedure 

Most of our sample consists of unresolved spectroscopic 
WD+dM binaries (see in Figure|5]for 16 representative spec- 
tra from our sample). To conduct a proper analysis of the 
individual stellar components, it is necessary to separate the 
individual spectra from the combined spectrum. To do this, 
we constructed an iterative method for separating the com- 
bined WD+dM spectrum into its i ndividual componen ts by 
fitting the app ropriate WD models (|Koester et al.l 120011) and 
dM templates dBocha nski et al. 20O7T) to each spectrum. The 
iKoester et all (120011) WD models consist of 372 DA and DB 
WD models for varying temperatures and log g. The DA mod- 
els vary from 6,000 to 100,000 K and each temperature bin 

9 http://vizier.u-strasbg.fr/ 

10 http://delia.phy. vanderbilt.edu/slowpokes?d=wdm 
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Figure 5. Examples of a typical SDSS WD+dM spectra reported in this 
paper. The fluxes have been normalized for visualization purposes. Note that 
large feature commonly found at around 5500A is an artifact caused by the 
SDSS spectrograph and is not a real feature (S06). 

has a separate model for log g values of 5, 5.5, 6, 6.5, 7, 7.5, 
8, 8.25, 8.5, and 9. The DB white dwarf models vary from 
10,000 to 40,000K and have logej values of 7, 7.5, 8, 8.5, 
and 9. We used half spectral type dM templates (M0, M0.5, 
Ml, M1.5...M8.5, M9; J. Bochanski, private communica- 
tion) that were derived from composites of high S/N spectra 
dBochanski et alj|2007h . 

The iterative process was comprised of five steps: 1) We 
began by fitting the combined spectrum (see Figure 5 for 
examples) with the WD models and subtracting the best-fit 
X 2 model from the combined spectrum; 2) we then fit a dM 
template to the WD subtracted spectrum to find the best-fit 
dM template; 3) then, we subtracted the best-fit dM tem- 
plate (from step 2) from the original combined spectrum and 
then found the best-fit WD to the dM subtracted spectrum; 
4) again, we found the best-fit dM template to the WD sub- 
tracted spectrum; and 5) we repeated this process until the 
best-fit WD and dM converged (did not change after repeated 
iterations). The process typically finished after only two or 
three iterations, however some fits were degenerate and alter- 
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Figure 6. Four examples of the results from our WD+dM spectral separation procedure as outlined in Section 3.1. Top panels show the original SDSS spectrum 
(black), with the recovered dM spectrum (orange; original SDSS spectrum subtracted by the best fit WD template), and the best fit dM template (red). On the 
bottom panel is again the original SDSS spectrum (black), with the recovered WD spectra (light blue; original SDSS spectrum subtracted by best fit dM template), 
and the best fit WD model (dark blue). 



nated between different values of a single parameter. For the 
few cases where solutions did not converge, we stopped the 
procedure after ten iterations and chose the model fit with the 
lowest overall x 2 value. We note that we do not distinguish 
between cold and hot WD solutions in our fitting process as it 
is not imperative for our magnetic activity analysis; this could 
be the cause of some of the solutions failing to converge in the 
fitting process. 

For the dM fit, we used a best-fit \ 2 technique over a wave- 
length range of 7000-9300A. We found that evenly weighting 
the entire wavelength regime of the dM provided the most ac- 
curate fits. However, the WD fitting process was more com- 
plicated; we placed higher weights on the hydrogen Balmer 
absorption lines and weighted them ~1000 times higher than 
the continuum to ensure that both the shape and the depth of 
the absorption lines were accurately fit, features important for 
determining \ogg. We focused on Ha, H/3, H<5, and Hy and 
weighted each predicted peak (corrected for RV of the WD, 
detailed below) in addition to 100 A region on either side of 
each peak to ensure the Balmer regions were fit well. Ad- 
ditionally, we weighted 10A wide portions between the ab- 
sorption lines ~1000 times higher to assist in properly fitting 
the overall shape of the WD continuum. This procedure re- 
sulted in accurate fits for both the Balmer lines as well as the 



continuum, which ensured an accurate subtraction of the WD 
from the combined SDSS spectrum and a clean dM spectrum. 
We also took into account possible hydrogen Balmer series 
emission from an active dM component that could affect the 
WD fit. We down-weighted the predicted locations of dM 
emission (corrected for RV) ~1000 times lower to minimize 
contamination. 

After each iteration, we used a cross-correlation technique 
to calculate the RVs of the WD and dM components. The 
cross-correlation technique compared the separated WD and 
dM to their corresponding rest wavelength best-fit template; 
the shift in the cross-correlation peak yielded the component 
RV. The measured RV was then used to shift the wavelength 
axis of the template and provide a better fit in subsequent iter- 
ations. For the cross-correlations, we used a wavelength range 
of 3900-5200A for the WD (this regime was chosen to include 
H/3, H 7 , and US) and 6600-9000A for the dM (this range in- 
cludes the strong Na iAA [8183.26, 8194.79A] doublet and 
a number of TiO molecular features). In an effort to deter- 
mine the uncertainty in the measured RV, we chose two bins 
over the wavelength ranges designated above for the WD and 
four bins for the dM component. For example, the WD had 
two bins 3900-4550 and 4550-5200A, while the dM had four 
bins: 6600-7200, 7200-7800, 7800-8400, and 8400-9000A. 
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Table 1 

WD+dM Parameters: Spectral classifications and radial velocities. 



SDSS ID 




Plate 


MJD 


Fiber 


dM 
SpT 


WD 
Type 


WDT e// 
(K) 1 


WD 

logg 1 


dM RV 

(km s ) 


WDRV 

n — 1 \ 

(km s ) 


SDSS J000109.424 


-255459.5 


2822 


54389 


334 


M2.5 


DA 


16000 


7.00 


-22.316 


208.306 


SDSS J000152. 104 


-000644.7 


387 


51791 


157 


MO 


DA 


40000 


7.50 


4.076 




SDSS J000250.65- 


045041.6 


2630 


54327 


439 


M4 


DB 


20000 


8.50 


-11.468 


-36.618 


SDSS JOO0356.93- 


050332.7 


2630 


54327 


173 


M4.5 


DA 


18000 


8.00 


-52.853 


73.650 


SDSS J000421.614 


-004341.5 


685 


52203 


561 


M5 


DB 


16000 


9.00 


-22.385 




SDSS J000447.784 


-291140.9 


2824 


54452 


600 


M0.5 


DA 


16000 


8.00 


-18.033 


-35.996 


SDSS J000453.944 


-265420.4 


2824 


54452 


78 


M4 


DA 


10000 


8.25 


-29.225 


174.521 


SDSS J000504.924 


-243409.7 


2822 


54389 


180 


M4.5 


DA 


11000 


8.25 


41.247 


29.640 


SDSS J00053 1.10- 


054343.2 


2624 


54380 


82 


M5 


DA 


12000 


7.75 


-23.974 


-46.497 


SDSS J000559.87- 


054416.1 


2624 


54380 


60 


M1.5 


DA 


30000 


8.00 


-9.603 


101.562 



We note that the values for WD T c ff and WD log g values are only estimates based on the best-fit models. The WD parameters reported 
n[2]aie calculated using methods described in Section 3.2 and are more reliable. 



For each ~600A bin, we shifted the wavelength regime by 
10A 30 times towards lower wavelengths (blue) and 30 times 
towards longer wavelengths (red) and calculated the RV us- 
ing the cross-correlation technique. The RV was taken to be 
the median value of all the measurements and we report the 
uncertainty as the median absolute deviation. The dM RVs 
are much more reliable because of the many narrow molec- 
ular features found in the dM spectrum, while the WD RVs 
have higher uncertainties due to their broad hydrogen Balmer 
absorption lines. 
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Figure 7. Distributions of various properties of the sample: a) the distri- 
bution of the dM spectral types of our sample (the black line is the entire 
sample, while the black dash-dotted line indicates the systems with high S/N 
that were used for the activity analyses); b) the distribution of the WD ef- 
fective temperature as derived by our model fitting procedure, which has a 
bimodal distribution at 10.000K and 19,000K; c) the distribution of derived 
WD log g values, which peaks around a value of 8.0; d) the distribution of 
the estimated binary projected physical separation in AU; e) the distribution 
of the distances (pc) calculated using the spectroscopic parallax relation; f) 
and distribution of calculated absolute distance above the Galactic plane (pc). 

Some examples of the results from the iterative fitting pro- 
cess can be found in Figure|6] The top panel of the plot shows 
the original SDSS spectrum in black, the derived dM flux in 
red (SDSS spectrum - best-fit WD model), and the best-fit 
dM template in orange. The bottom panel is similar: the orig- 
inal SDSS spectrum is shown in black, the derived WD flux 
in light blue (SDSS spectrum - best-fit dM template), and the 
best-fit WD model in dark blue. Our procedure returned dM 
spectral type, WD effective temperature (T er r), WD surface 
gravity (logg), and WD and dM RV . Histograms of these 
parameters can be found in Figure [7] and the values can be 



found in Table Q] We note that the WD parameters calculated 
through the above analysis are only estimates. In the follow- 
ing section (Section 3.2) we use the resulting WD spectrum 
(from the WD+dM separation procedure described above) to 
perform a more detailed analysis to calculate WD parameters . 
When available, calculated WD parameters reported from the 
analysis discussed in Section 3.2 and found Table [2] should be 
used. 

3.2. WD cooling analysis 

We derived the atmospheric parameters of the DA WDs 
(T e ff and \ogg) by performing a fit of the observed Balmer 
lines to the WD models of D. Koester (private commu- 
nication). The Balmer lines in such WD models were 
calculated with the modifie d Stark broadening profiles of 
iTremblav & Bergeronl (120091) . kindly made available by the 
authors. For the line fitting we used the code fitsb2, 
which follows a proc edure based on \ 2 minimization 
(iNapiwotzki et al]l2004l) . We considered the lines from H/3 
to He for the fitting procedure when S/N<30 and from H/3 
to H£ when S/N>30. We did not include Ha to minimize 
any effects from the M dwarf spectral subtraction technique 
(Section 3.1). As initial values, we considered the photomet- 
ric temperature and a log g = 8. The photometric T e g were 
obtained by fitting the S PSS synthetic gri magnitude s to WD 
synthetic photometry of iHolberg & Bergeronl (120061) . Since 
the SDSS gri magnitudes are likely contaminated by the M 
dwarf component we derived synthetic gri magnitudes from 
the spectrum of the WD component separated by the proce- 
dure outlined in Section 3.1. The entire WD fi t ting pr oce- 
dure is descr i bed in more detail in IGarces et al.l (1201 ll) and 
iCatalan et alj d2008bl) . By using the photometric T c ff as ini- 
tial value in the Balmer line fitting procedure we avoided 
possible degeneraci es that result from the cool/hot solution 
(IGarces et al.l 1201 lh . In Figure [8] we show the fits to the 
Balmer line pr ofiles for 15 DA WDs in our sample. 

According to lBergeron et al.l(|1992l) . the atmospheres of DA 
WDs below 12000K could be enriched in He while preserv- 
ing their DA spectral type. This He is thought to be brought to 
the surface as a consequence of the development of a H con- 
vection zone. Depending on the efficiency of convection, the 
star could still show Balmer lines, instead of being converted 
into a non-DA white dwarf. The increase in He mimics the 
behavior of a high gravity, thus, the log g obtained from the 
spectroscopic fitting can not be trusted. For this reason when 
T G ff <12000K we prefer to consider the photometric T c s and 
the typical value for the surface gravity logy = 8.0. 
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Table 2 

WD+dM Parameters II: WD Cooling analysis parameters 



SDSS ID 



WDT cff (ffWDT„ f[ ) 
(K) 



M WD (o"M WD ) 
(M ) 



T c ool ( CT T cool ) 
(Myr) 



SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 



J000152.10- 
J000356.93- 
J000453.94- 
J000504.92- 
J000531.10- 
J000559.87- 
J000651.91- 
J000711.41- 
J000829.92- 
J 001247.18 



-000644.7 
-050332.7 
-265420.4 
-243409.7 
-054343.2 
-054416.1 
-284647.1 
-283848.6 
-273340.5 
1-001048.7 



7598.0(166.0) 

17617.0 (453.0) 
9503.0 (285.0) 

12916.0 (594.0) 
13753.0(1505.0) 

29981.0 (523.0) 
14558.0 (1609.0) 
13401.0(1854.0) 

17068.0(801.0) 
18843.0(2922.0) 



0.594 
0.661 
0.599 
0.826 
0.789 
0.767 
0.511 
0.732 
0.451 
0.993 



(0.0006) 
(0.0464) 
(0.0080) 
(0.0627) 
(0.1019) 
(0.0532) 
(0.1057) 
(0.1254) 
(0.0786) 
(0.2610) 



1150.300(47.500) 
132.119(14.636) 
651.520(37.821) 
490.303 (78.049) 

373.370(105.280) 
16.788 (3.645) 
160.609 (51.463) 

350.238 (119.961) 
77.320(16.291) 

289.300 (201.500) 



2.0 



CD '-5 

> 1.0 



2.5 
2.0 









T e „= 28870 + 430 K 
log g- 8.07 ± 0.09 


T 8 „= 24780 4 770 K 
log g- 7.38 4 0.1 1 


T e „= 15890 4 420 K 
log g= 7.84 i 0.10 








T 8 „= 27900 4 470 K 
log g= 7.79 ± 0.08 


T.„= 20500 -+ 510K 
log g= 7.78 ± 0.09 


T e „= 12370 4 370 K 
log g= 8.02 ± 0,12 








T„„= 23410 ± 760K 
log g= 7.43 ± 0. 1 1 


T,„ = 19640 ± 480 K 
log g= 7.92 ± 0.09 


T.„= 34630 ± 510K 
log g= 7.89 ± 0.09 








T„„= 30240 ± 410 K 
log g= 7.39 ± 0.09 


T„„= 21790 ± 520 K 
log g = 7.85 ± 0.09 


T„.,= 16320 ± 400 K 
log g= 7.97 ± 0.09 








T„„= 22600 ± 580 K 
log g= 7.88 ± 0.09 


T„„= 29570 ± 390 K 
log g= 7.78 ± 0.08 


T„= 12760 ± 380 K 
log g= 8.00 ± 0.02 
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Figure 8. Fits to the individual Balmer lines for the WDs for 15 separated 
WD in our sample. The line profiles correspond to H/3 (bottom) up to He or 
(top). We have applied vertical shifts for clarity. The solid lines are the 
observed spectra and the red lines are the models that best fits the profile. 

Once we derived the T c g and log g for each WD, its mass 
(Mwd) and cooling time (t coo \) were obtained using ap- 
propriate co oling sequences. W e have adopted the cool- 
ing tracks of ISalaris et al.l (120001) . which consider a carbon- 
oxygen (C/O) core WD (with a higher abundance of O at the 
central core) and a thick hydrogen envelope on top of a helium 
buffer, g(H) = M H /M = 10~ 4 and g(He) = M He /M = 
10~ 2 . The resulting values for Mwd and t coo \ are listed in Ta- 
ble |2] The total age of a WD is comprised of its cooling time 
plus the lifetime of its progenitor star during the pre-white 
dwarf stage. If a white dwarf is isolated, or in a wide binary 
(no interaction between members), we can calculate the mass 
of the progenitor usi ng an initial-final mass relationship (e.g., 
iCatalan eta l. 2008a), and considering stellar tracks, estimate 
the progenitor lifetime. However, in this study we are consid- 
ering white dwarfs in close binaries that may have interacted 
in the past, transferring mass and altering the estimated pro- 
genitor properties. However, the cooling time of the WD can 
be used to constrain its total age by setting a lower bound. 
This lower limit is particularly important when the WD com- 



panion is an dM, si nce determining t he ages of dMs is partic- 
ularly challenging (ISoderblomll20 1 01) . We derived T e g, log <j, 
Mwd, and i coo i and their corresponding uncertainties for 904 
WDs in our WD+dM sample. The values are reported in Ta- 
ble m 

3.3. WD+dM parameters 

In addition to the values presented in Tables Q] and [2] we 
derived additional parameters that are given in Table [3] In 
particular we derived the system RV through space (km s" 1 ), 
the height above the Galactic plane (pc), approximate distance 
to the system (pc), approximate binary separation (AU) and 
orbital period (days) assuming circular, Keplarian orbits and 
an edge-on inclination, activity state (as traced by Ha), and 
magnetic activity strength (using Lna/Lboi)- In this section 
we discuss in detail on how we calculated each of these pa- 
rameters. 

3.3.1. System radial velocity 

In our sample, 1714 dMs and 1674 WDs had reliable RV 
measurements, resulting in 1627 WD+dM pairs with RV 
measurements for both components. In addition, we derived 
WD masses for 904 pairs as discussed in Section 3.2; for 
the rest of the WDs, we assigned a mass of 0.6 M B (the 
avera ge mass of 90% of WDs in SDSS DR4; iKepler et alj 
|2003). We inferred masses of the dMs from their spectral 
types dReid & Hawlevlr2005b . Then using center-of-mass ar- 
guments in circular orbits, MdM^'dM = Mwd^wd, along with 
the measured relative RVs, we estimated the system veloc- 
ity of the pair. The measured RVs are relative velocities ex- 
pressed as w wd = u w( j + w sys , where w wd is the measured RV 
of the WD, w sys is the system velocity, and u w d is the abso- 
lute RV of the WD. Similarly, we measured w dM for the dM: 

v dM = ^dM+^sys- Combining these three equations, we solved 
for the system velocity: 

Absolute WD and dM RVs were calculated using the mea- 
sured relative RVs and the derived system velocities. Absolute 
WD and dM RVs, system velocities, and their corresponding 
uncertainties are reported in Table [3] uncertainties were cal- 
culated using standard error propagation. We note that our 
reported system velocities are lower limits due to the assump- 
tion of circular, edge-on, and location in the orbit. In reality, 
most of our measured RVs have lower values than the true 
stellar orbital motions. Thus, the system velocities are used 
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for statistical relative comparisons rather than exact measure- 
ments for individual pairs. 

3.3.2. Separations & Periods 

Without multi-epoch RV measurements, system parameters 
could not be determined via full orbit solutions. Instead, we 
assumed the single-epoch absolute RV of each component (or 
V — r sin i) to be its orbital velocity and calculated the pro- 
jected linear separation for the binary assuming a Keplerian 
circular orbit and edge-on inclination. The total separation of 
a system was found by summing the individual component 
separations from the center-of-mass of the system, which was 
calculated by balancing centripetal and gravitational force. 
Considering our calculated RVs are lower limits, this sug- 
gests that our measured separations are likely upper limits. 
When multiple RV measurements were available (see Sec- 
tion 3.4), we constrained the limits on RVs and separations 
using the set of WD and dM measured RVs that had the max- 
imum difference between them; this is where the RVs will be 
most representative of the orbital velocity. While the projec- 
tion effects can be statistically corrected for ensembles, it is 
a steep function of o rbital separation and total system mass 
(iDupuy & Liull20i"il) . Applying the correction for individual 
systems introduces large uncertainties and does not affect the 
relative scalings; hence, we chose to not use the corrections. 
The estimated projected linear separations (AU) and orbital 
periods (days) are reported in Table [3] 

3.3.3. Distances 

We used a spectroscopic par allax relation that was derived 
from data in IBochans ki (2008) to estimate the absolute mag- 
nitudes of the dM as a function of spectral type and thus calcu- 
late the approximate distance to the dM. We estimated that the 
distances are accurate to 20%. Photometric parall ax relations 
for dMs are more accurate measures of distance (IWest et al.l 
12001 iDavenport et all l2006t IBochanski et al.l 1201 a 1201 11) . 
however, due to the photometric contamination by the close 
WD companion (in the r-band), the resulting distances would 
have significant errors. Calculated distances are reported in 
Table [3] and a distribution of distances can be found in Fig- 
ure |7] 

3.3.4. Vertical Distance from the Galactic plane 

Using the derived distances from Section 3.3.3, the Galac- 
tic latitude and longitude, an d by assumin g the Sun is 15 
pc above the Gala ctic plane dCohe 



iBinnev et al 



assuming 
n 1995; 



_ Ng et al.lfl997l: 
119971) and at a distance of 8.5 kpc from Galactic 



center (e.g., iMatsunaga et al.ll2009l) . we calculated the abso 
lute height above the Galactic plane for each pair using the 
following relation: 



Gal 



d x sin(fe — atan 



*cos(/)) + z Q , (9) 



where Za a i is the height above the Galactic plane, I and b are 
galactic latitude and longitude, r© is the Sun's distance from 
Galactic center (pc), z is the Sun's height above the Galactic 
plane, and d is the distance measured from Section 3.3.3. We 
used the height above the Galactic plane as a statistical proxy 
for the age of the binary system; younger stars are formed in 
the plane of the Galaxy and ove r time are dynam ically heated 
and perturbed to heigher orbits (IWest et al.ll2006t W08). Dis- 
tributions of calculated Galactic heights are reported in Fig- 
ure|7]and the values can be found in Table|3] 
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Figure 9. Examples active dMs (left panel) and inactive dMs (right panel) 
from our WD+dM sample. The red dot-dashed line is the location of the 
radial velocity corrected Ha line as measured by our WD+dM separation 
procedure outlined in Section 3.1. The Ha EW and S/N are printed on each 
panel. 

3.3.5. Magn etic Activity 

We use the Ha emission line as an indicator for magnetic 
activity in the dMs. In order to calculate the magnetic activ- 
ity strength, the line flux of Ha was measured using a trape- 
zoid al integration tec hnique over 12 A centered on the Ha 
line (I West et al.ll2004l) . The continuum was taken as the aver- 
aged flux 5A on either side of the line region, the noise was 
the standard deviation of the flux over the continuum range. 
We calculated t he Lna/^bo\ of each d M using the meth- 
ods de scribe d bvlWalkowicz et al.l (120041) and lWest & Hawlevl 
(120081) . In IWalkowicz et al.l (120041) . the authors derived a 
spectral type dependent scaling factor, \ (the ratio of contin- 
uum flux near Ha over the bolometric flux), from a sample of 
nearby stars and 2MASS standards. Multiplying the spectral 
type dependent x values by the measured Ha equivalent width 
(EW) resulted in Lna/Lboi independent of the dM distance. 
Lna/Lboi is a measure of the flux output driven by magnetic 
activitiy via Ha emission over the total bolometric luminosity 
of the star; this quantity is independent of the stellar contin- 
uum (unlike EW). Magnetic activity strength as measured by 
Lna/Lboi is reported in Table [3] 

The magnetic activity state of the dMs was determined by 
inspecting the Ha line of each dM (WD subtracted) spectrum 
by-eye and separating them into three main categories; ac- 
tive, inactive, or unable to classify due to low S/N spectra. To 
ensure the high quality of our classifications we only consid- 
ered objects to be active or inactive if the S/N near Ha was 
greater than three. Of our WD+dM sample, we classified 741 
as active and 496 as inactive - 526 were unclassified. Four of 
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Table 3 

WD+dM Parameters III: Additional parameters 



SDSS ID 




dMRV 1 


WDRV 2 


System RV 3 


z 


d 


s 4 


Period 


Active 


LHn/Lbol 






(km s" 1 ) 


(kms" 1 ) 


(kms" 1 ) 


(pc) 5 


(pc) 6 


(AU) 7 


(days) 






SDSS J000109.424 


-255459.5 


-22.3 (2.6) 


208.3 (208.9) 


116.1 (125.7) 


-223.6 


410.3 


0.08 


8.5 


no 




SDSS J000152.104 


-000644.7 








-894.6 


1047.7 






yes 


1.04e-04 


SDSS J000250.65- 


045041.6 


-11.5 (8.4) 


-36.6 (26.5) 


-28.2(17.7) 


-275.0 


320.3 


8.54 


9626.5 


yes 


1.63e-04 


SDSS J000356.93- 


050332.7 


-52.9 (6.3) 


73.6 (16.0) 


33.7 (14.5) 


-121.8 


150.7 


0.40 


93.8 


yes 


1.82e-04 


SDSS J000421.6H 


-004341.5 








-308.2 


373.5 






no 




SDSS J000447.784 


-291140.9 


-18.0 (5.9) 


-36.0 (25.9) 


-27.0(13.1) 


-2070.1 


3876.1 


13.23 


16066.7 


no 




SDSS J000453.944 


-265420.4 


-29.2(11.4) 


174.5 (48.2) 


106.6 (33.9) 


-192.8 


364.3 


0.13 


18.0 


yes 


1.24e-04 


SDSS J000504.924 


-243409.7 


41.2(4.1) 


29.6(14.5) 


32.7 (8.3) 


-75.7 


150.5 


80.3 


248091.09 


yes 


2.76e-04 


SDSS J00053 1.10- 


054343.2 


-24.0 (3.7) 


-46.5 (85.9) 


-40.8 (64.4) 


-33.6 


53.2 


21.4 


36401.82 


yes 


9.88e-05 


SDSS J000559.87- 


054416.1 


-9.6 (2.3) 


101.6(35.7) 


60.3 (23.5) 


-471.2 


532.1 


0.5 


110.13 


yes 


8.82e-05 


1 The absolute dM radial velocity corrected for the system velocity. The 


corresponding dM RV uncertainties are repoited in the 


parenthesis. 







The absolute WD radial velocity corrected for the system velocity. The corresponding WD RV uncertainties are reported in the parenthesis. 
3 The system velocity as defined in Section 3.3.1. The system RV uncertainties calculated from the dM and RV uncertainties using standard error propogation techniques 
are reported in the parenthesis. 

^ Projected linear separation as defined in Section 3.3.2. Recall that the separations presented in this paper are lower limits. 
The absolute height above the Galactic plan calculated using the Equation 9 

6 Distances were calculated using a spectroscopic parallax relation derived from data found in Bochanski 1 2008). 

7 The separations were calculated using the dM and WD RVs after being corrected by the system velocities. We assumed that the measured absolute RVs of each 
component is the orbital velocity as well as assuming edge-on circular Keplarian orbits. 



our objects also have very broad Ha emission f eatures, possi- 
bly due to accretion (e.g.. lHartigan et al. 19951; we classified 
these separately and excluded them from our magnetic activ- 
ity analysis. A few typical examples of active and inactive 
stars can be found in Figure|9] The activity states are reported 
in Table[5]as "yes" for active stars, "no" for inactive stars, and 
"N/A" for unclassified stars. 



3.4. Multiple exposure RV analysis 

All SDSS spectra are composite spectra made up of mul- 
tiple exposures (at least three but as many as 15). The typi- 
cal exposure times were roughly 900-1200 seconds and were 
often made in succession. These individual exposures have 
been used in previous s tudies to investigate M dwarf flare 



rates ([Hilton et al 



dKruse et alj|201 



20101) and time variable spectral features 



Bellet alj|2012ri . We used the multiple ex 



posure spectra in an effort to investigate time variability in the 
derived radial velocity measurements in our close WD4-dM 
binary pairs. Many of our pairs are very closely separated, and 
thus may have periods on the order of a few hours. The single- 
epoch radial velocities were calculated in a similar manner as 
discussed in Section 3.1, except that the best-fit template for 
each spectrum was already determined. We also calculated 
LHa/Lboi for each individual exposure as described in Sec- 
tion 3.3.5. The number of exposures, the median change in 
radial velocity (ARV), the median error in the change of ra- 
dial velocity (ctarvX an d the time in between exposures of 
median RV change (hours) are reported in Table |4] 

Being able to measure multiple RVs during the orbit al- 
lowed us to put better constraints on the separation of the 
WD+dM pairs as there was a greater chance of catching the 
components at their maximum separation in one of the expo- 
sures (see Section 3.3.1). Very close pairs are apparent by the 
rapidly changing radial velocities over the baseline of the con- 
secutive observations. In addition, these very close pairs may 
be interesting for followup observations and may be used to 
investigate the variability of magnetic activity as a function of 
orbital motion. 



4. RESULTS 

4.1. Activity fraction vs. Galactic height 

Figure [10] shows the activity fraction in close WD+dM 
pairs (filled circles) as a function of spectral type compared to 
the activity fractions in the dM field sample (open diamonds, 
Wll). The uncertainties shown in Figure [10] are calculated 
from the binomial distribution. We see a noticeably higher ac- 
tivity fraction in the close WD+dM pairs than in the isolated 
field dMs across all spectral types. However, towards later 
spectral types we see the disparity in the activity fractions de- 
crease between the WD+dM and isolated dM populations. 
The convergence of the two populations at late spectral types 
seems to indicate that having a close binary companion may 
not affect activity in late-type M dwarfs in the same manner 
that it does for the earlier-type dMs. 
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Figure 10. Fraction of active stars vs. spectral type for our WD4-dM sample 
(solid circles) and isolated field M dwarfs (open diamonds, Wll). The error 
bars were computed from the binomial distribution. The M6 and M7 bins 
were grouped together to help increase the counting statistics for the later 
types. The fraction of active WD4-dM is significantly higher than that of 
field M dwarfs. Towards later spectral types we begin to see a decrease in the 
difference in activity fractions. 
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Figure 11. Fraction of active stars as a function of absolute Galactic height (pc) for each spectral type (proxy for age, W08). The filled circles are bins of 
WD+dM (a bin is required to have two or more stars) and the open triangles are bins of field dMs (Wll). In an effort to increase the binomial statistics in the 
least represented bins, we combined the M0 and Ml spectral types as well as the M6 and M7 spectral types. There are three main trends: 1 ) M dwarfs in WD+dM 
pairs are more active (up until M5) than their single M dwarf counterparts; 2) whatever the mechanism driving the dynamo in fully convective late-type dMs, it 
does not seem to be affected by a close binary companion (or possibly rotation); and 3) we see elevated activity fractions close to the plane and the presence of 
activity decreasing farther from the plane (age), indicating lifetime in early-type M dwarfs with WD companions is most likely finite but lasts longer than field 
dwarfs. In late-type M dwarfs with WD companions the lifetimes are comparable to the late-type field population activity lifetimes. 



Some of the morphology of Figure QJjis due to age effects 
(early-type M dwarfs have shorter active lifetimes than their 
late-type counterparts; W08) and the location of stars in the 
Galaxy (which is age dependent). Figure [TTJ puts our sample 
in the proper Galactic context by showing the activity fraction 
as a function of absolute Galactic height (a proxy for age) for 
each spectral type. The WD+dM (filled circles) sample is 
binned by Galactic height and we chose our bin sizes to opti- 
mize the number of stars per bin; the typical bin (filled circles) 
has approximately 10 stars (with a minimum of 2 stars). We 
used the same bin sizes for the field M dwarfs (open triangles, 
Wll). Due to the relatively small numbers of some spectral 
types, we combined M0 with Ml and M6 with M7 for our 
analysis. There are no M8 or M9 dwarfs with reliable activity 
parameters that were identified in our analysis. 

FigurefTTIshows that in comparison to isolated dMs (Wl 1), 
the activity fraction is noticeably higher in almost all Galac- 
tic height bins until a spectral type of M5, confirming the 
higher activity fraction trend seen in Figure QJj] Also, at a 
spectral type of M5, the disparity in activity fractions between 
field dM and close WD+dM pairs decreases and continues to 
M6+M7 bins. Due to the high uncertainties the trend is not 
as clear in the latest type bins. Coincidently, this change oc- 
curs roughly at the spectral types with interiors that are fully 
convective. In addition, the activity fraction as a function of 
Galactic height can be interpreted as evidence for activity life- 
time. By comparing the active fractions of dMs to simple 
ID dynamical simulations, W08 showed that early-type dMs 
have short activity lifetimes (~0. 6-2.0 Gyr for M0-M3, re- 



spectively) and that late-type dMs have longer activity life- 
times (^7.0-8.0 Gyr); there is a sharp increase of activity life- 
times at a spectral type of M4 (~4.5 Gyr), the location of the 
onset of full convection in dMs. In Figure QT| the low ac- 
tivity fraction in the isolated early-type dMs (Wll) can be 
explained by the short activity lifetimes described by W08; 
the activity has already had enough time to "turn off" near the 
plane. However, the isolated late-type dMs (Wll) have high 
activity fractions close to the plane and continue to be active 
farther from the plane, indicating these stars remain active for 
longer. When examining the close WD+dM population, Fig- 
urefTTIshows a larger activity fraction (when compared to the 
field population) of early-type WD+dMs close to the plane 
and an elevated activity fraction decreasing steadily farther 
from the plane, whereas almost no early-type field dMs were 
active farther from the plane. These results provide evidence 
for extended yet finite activity lifetimes for early-spectral type 
WD+dMs (M0-M3). For the late-type dMs, FigurefTTIshows 
that the close WD+dM activity fractions are similar to the 
late-type field dM population close to the plane and continue 
to be comparable at higher Galactic heights. This similarity 
is interpreted as the late-type close WD+dMs and late-type 
field dMs having the same or similar activity lifetimes. Coin- 
cidentally, M5 is just beyond the spectral type where the onset 
of full convection occurs in the interior of dMs. 

4.2. Activity strength vs. Galactic height 

Figures [TOl and [TTI indicate that WD+dM with early-type 
dMs are more likely to be active than their field counterparts 
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Figure 12. Activity strength (as measured by Lna/Lboi) vs. dM spec- 
tral type. The data are binned by spectral type and the median value of 
Llia/Lbol i s plotted. The error bars represent the 25th and 75th quartiles 
in that bin. Red filled circles represent binned WD+dM and blue open trian- 
gles represent isolated dM from Wll. We see a decrease in activity strength 
for field dwarfs at late types as well as a decrease in strength for paired 
dwarfs. The difference in activity strengths between the two populations in 
the M6+M7 bin indicates that having a WD companion effe cts the activity 
strength but not the activity fraction as was seen in Figure 1101 



while late-type dMs in close WD+dM pairs show little differ- 
ence in their activity state compared to late-type field dMs . 
Figure [reexamines the median activity strength (as measured 
by L Hct /Lboi) of the M dwarfs in WD+dM (filled red circles) 
in comparison to the median activity strength in field dwarfs 
(open blue triangles, Wll). The error bars represent the 25th 
and 75th quartiles of Lna/Lboi in each bin. The activity 
strength in WD+dM pairs follows the same trend that we see 
in the field dM sample (Wll); the activity strength is constant 
in early spectral types until a spectral type of M5, when the 
activity strength decreases. The median activity strength in 
the WD+dM pairs is slightly larger than the median activity 
strength for the field dM population in every spectral type bin. 
Similarly, the amount by which Ljja/Lboi in WD+dM pairs 
is larger than field dMs remains constant (log(LHa/Lboi)~ 
+0.1 dex) until a spectral type of M6+M7, where the differ- 
ence in Lna/Lboi between WD+dM pairs and field dMs be- 
comes noticeably larger (log(LH a /Lboi)~ +0.3 dex). This 
suggests that there is an increase in activity strength for later 
spectral types with close companions. Unfortunately, we 
lack large numbers of late-type dMs (M8 and later) to ex- 
plore this trend further. To determine how similar the two 
(WD+dM and isolated dMs) populations are to one another, 
we performed a tw o-sided Kolmogorov-Smirnov test (KS test, 
iPress et alll 1992b to determine whether the Lna/Lboi for all 
the spectral type bins between the WD+dM and isolated dM 
populations can be drawn from the same parent population. 
We found an 87.1% chance that the two populations are sig- 
nificantly different. Figure Q~2] demonstrates that the activity 
strengths, along with the activity fractions, are consistently 
stronger for dMs with close WD companions. 

Analogous to Figure[JTJ Figure[l3]investigates how the me- 
dian activity strength (traced by Lna/Lboi) behaves as a func- 
tion of Galactic height and spectra l typ e. We retain the same 
binning and formatting for Figure [13] that was used for Fig- 
ure QTJ Figure Q~3] shows that the activity strength of the dMs 
in WD+dM pairs remains consistently higher than the activ- 
ity strength of field dMs with increasing Galactic height for 
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Figure 13. Activity strength (as measured by Lna/Lboi) as a function of 
Galactic Height for different dM spectral types. The filled red circles are 
the median Ljia/Lbol value of the bins of WD+dM, a bin is required to 
have two or more stars. Isolated M dwarfs from Wll are represented in 
blue open triangles. The error bars represent the 25th and 75 th quartiles of 
Lija/Lbol f° r eac h Dm - The WD+dM exhibit higher activity strengths than 
the field population across all spectral types. From spectral types from MS- 
MS there is a slight decrease in the dM field populations Ljia/Lbol with 
increasing Galactic height, whereas, there appears to be no such decrease 
with the WD+dM population. Also, at later-types (M6+M7) there is a larger 
difference in the Lna/Lboi values between the two populations than seen in 
earlier spectral types. A two-sided KS test found high significance that these 
two (M6+M7) populations are drawn from a separate parent population. 

spectral types M0-M4. We calculated the KS test for M0+M1, 
M2, M3, and M4 spectral type bins and found that there is 
a 98.3, 66.2, 99.6, and 97.8% chance, respectively, that the 
binary and field populations are drawn from different parent 
populations. At spectral types > M5, the activity strength 
in field dwarfs decreases with increasing Galactic height (or 
age), a trend that is not as obvious in the paired dMs. Using 
the KS test we report a 99.8% likelihood that the two M5 pop- 
ulations are drawn from a different parent population. Simi- 
lar to M5, we see evidence that dMs in WD+dM with spec- 
tral types M6+M7 had larger activity strengths than the field 
dwarf counterparts. With only three bins, we were unable to 
verify the decreasing activity strength with increasing Galac- 
tic height as seen in the M5 bins. Since we only have three 
bins, we performed a KS test for the stars within each of the 
three bins. In order of increasing Galactic height, we found 
98.9, 99.5, and 78.7% significance that the two populations 
are drawn from different populations. Given the paucity of 
stars of our sample in the later spectral types, we can only 
speculate that the activity strength appears to remain strong 
in the later spectral types and also possibly with increasing 
height (or age). Figure Qj] suggests that there may be a shift 
in the behavior of activity strength (as a function of spectral 
type and Galactic height) at spectral type M5, coinciding with 
a change in the structure of dM interiors. 

4.3. Activity & Activity Strength vs. Binary Separation 

In Figure [14] we plot the fraction of active stars as a func- 
tion of the dM spectral type for three separate separation bins 
using our projected linear separations (Section 3.3.2): < than 
0.1 AU, 0.1-1 AU, and 1-100 AU. Stars with separations < 0.1 
AU (purple solid line, diamonds) tend to be more active than 
stars with separations between 0. 1-1 AU (blue dotted line, tri- 
angles), which are in turn more active than pairs with sep- 
arations between 1-100 AU (red dashed line, squares). We 



14 




M2 M3 M4 
dM Spectral Type 



M5 M6+M7 



Figure 14. Activity fraction for WD+dM versus spectral type grouped by 
projected, linear separation: < 0. 1 AU (puple solid line, diamonds), from 
0.1-1 AU (blue dotted line, triangles), 1-100 AU (red dashed line, squares), 
and the Wll field dM population (black dash-dotted, solid circles). In gen- 
eral, closer separated pairs tend to be more active than wider pairs, and even 
the widest pairs have elevated activity levels compared to the field popula- 
tion. The narrowing of the activity fractions in the three separation bins to- 
wards later spectral types is significant and may provide evidence that close 
companions no longer effect the activity state of late-type dMs. 

also include the Wll (black dash-dotted line, solid circles) 
isolated dM sample to compare activity behavior between the 
most widely separated pairs and the field sample. Figure [14] 
demonstrates that at early spectral types, the widest pairs are 
unlikely to be active and have similar activity fractions as the 
field dM population. Towards later spectral types, the wide 
pairs appear to be more active than isolated stars, however, 
this result is less robust due to the increasing uncertainty in 
the late spectral type bins for the widest pairs. Figure [Pflalso 
shows that the difference in activity fractions between the sep- 
aration populations decreases with increasing spectral type. 
We see a possible convergence of activity fraction in the later 
spectral type bins of M6+M7. Once again, this convergence 
appears to begin at the M3/M4 transition to fully convective 
interiors, which may indicate that a different mechanism is 
responsible for generating the magnetic fields in late-type M 
dwarfs. 

Figure Q3] shows the log(LH a /Lboi) vs. log(binary separa- 
tion) [AU] for all active pairs broken down by spectral type. 
For each spectral type bin we performed a linear least-squares 
fit to the data; the slope of each fit can be found on the corre- 
sponding panel in Figure [15] We also performed a goodness- 
of-fit test (F Test) comparing the linear least-squares fit to a fit 
using a constant line through the data. The significance val- 
ues of the linear least-squares fits (as compared to a flat hor- 
izontal line) were 53.2, 11.3, 23.8, 1.9, 10.1, 53.5, 54.2, and 
31.0% for spectral types M0, Ml, M2, M3, M4, M5, M6, and 
M7, respectively. Our results suggest that there is no conclu- 
sive evidence that the binary companion increases the mag- 
netic strength of the dM in a way that correlates with binary 
separation. We found no distinct change in log(LHa/Lboi) 
vs. log(binary separation) across the M3/M4 fully convective 
boundary. From Figures [l4]and[T5] we see that the proximity 
between the WD+dM pairs does not appear to strongly affect 
the strength of the dM chromospheric activity. 

4.4. WD Cooling track ages 

Using the analysis described in Section 3.2 we were able to 
derive physical properties for 904 of the WDs, including T e s, 
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Figure 15. log(LH Q /Lt,oi) vs. log(binary separation) for active dMs in 
close WD+dM pairs for each spectral type. The solid line shows the linear 
least-squares fit for each distribution. All of the spectral types exhibit a nega- 
tive slope (with the exception of the M7, an inconsistency we attribute to the 
low numbers of pairs), implying higher magnetic activity strengths in the dM 
for smaller binary separations. However, as discussed in Section 4.3, the lin- 
ear least-squares fits are not statisticaly significant and the data are consistent 
with a horizontal line. 

mass, and cooling age (and their associated uncertainties). 
The WD masses were used to calculate the system velocities 
as well as the separations and periods of the pairs. We used 
the WD cooling ages as an alternative proxy for the age of the 
WD+dM systems. The WD cooling ages are accurate to a few 
100 Myr for WDs with T eff < 25000K and accurate to ~20% 
for WDs > 25000K; the WD cooling ages do not take into 
account the progenitor lifetimes. Therefore, the cooling ages 
are used in this study as a lower limit of the age and as a rough 
guide to distinguish between old and young systems. The pro- 
genitor ages cannot be determined for these systems as many 
of these pairs have close enough separations that mass transfer 
may have occurred; a slight change in the WD mass will have 
a significant effect on the es timated progenitor m ass (and its 
main sequence lifetime; e.g JCataran et alJ200 8b). Figure[T6l 
shows the WD cooling ages (Gyr) as a function of dM spec- 
tral type for the the active pairs (top panel) and the inactive 
pairs (bottom panel). WDs with log g errors > 0.4 and S/N 
< 5 were removed because of the large uncertainties in the 
calculated WD cooling ages. Figure [16] shows that there is 
dearth of early-type WD+dM pairs with active dMs in com- 
parison to the inactive dMs. Past spectral types of M4, both 
active and inactive pairs have a similar spread in WD cooling 
ages. We performed KS tests for spectral type bins M0+M1, 
M2, M3, M4, M5, and M6+7 and found that there is a 2.5, 
1.7, 2.5, 82.3, 46.9, and 100.0% possibility, respectively, that 
the active and inactive populations are drawn from the same 
parent population. From the results of the KS test, we con- 
clude that early-type, active dMs in WD+dM pairs tend to 
have younger WD companions than the inactive sample. We 
see less of a difference in the active and inactive populations 
in the late-type dMs. In the S06 study (in which they also used 
WD cooling ages), WD+dM pairs with active dMs in well 
sampled spectral type bins tended to be younger on average 
than WD+dM pairs with inactive dMs. Our results argue that 
only the early-type (< M3) active dMs appear to be younger 
on average. 

We note that there were three active, early-type WD+dM 
pairs that exhibited old WD cooling ages, which was incon- 
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Figure 16. M dwarf spectral type vs. WD cooling age (Gyr) as a function 
of activity state. The top panel is the WD+dM pairs with active dMs and 
the bottom panel is the WD+dM pairs with inactive dMs. There is a dearth 
of old active early-type dMs in WD+dM pairs indicating that dM activity 
lifetimes are short-lived as is seen in the field dM population (W08). By 
visual inspection there is no difference in the WD cooling ages spanned by 
the active and inactive late-type dMs. A KS test for the spectral type bins 
M0+M1, M2, M3, M4, M5, and M6+7 found that there is a 2.5, 1.7, 2.5, 
82.3, 46.9, and 100.0% possibility, respectively, that the active and inactive 
populations are drawn from the same parent population. Thus, for early-type 
WD+dM pairs, the active population is distinct from the inactive population, 
but there is no such difference for the late-type WD+dM populations. S06 
conducted a similar study and also found a dearth of active early-type dMs in 
WD+dM pairs. 

gruent with the rest of the sample. We suspect that the two 
pairs in the M0/M1 and M2 spectral type bins were active 
due their very small separations (~ 0.05AU; see Figure IT4Y 
The WD+dM pair in the M3 spectral type bin did not have 
as small of a separation (~ 0.17 AU), however, it is pos- 
sible that an inactive dM was observed during a flare event 
(iKowalski et alj|2009h . 

In Figure [17] we plot the WD cooling age against the abso- 
lute Galactic height (Figure [T7J left panel) and distance (Fig- 
ure [TTJ right panel). Both the WD cooling ages and Galactic 
heights are proxies for the age of the binary system and we 
expected to see a linear relationship in Figure [17] (left panel), 
albeit with some scatter. However, instead of the expected lin- 
ear relationship we identified a potential selection bias; there 
is an overdensity of systems with small WD cooling ages that 
span a large range of Galactic heights. Figure [l7]suggests that 
our sample is biased towards systems that are younger; young 
WDs are intrinsically brighter and are easier to observe farther 
away. Due to the SDSS footprint (which focuses on the North- 
ern Galactic cap), distance and absolute Galactic height are 
roughly equivalent. The similarity between Galactic height 
and distance can be seen in similar trends between the left 
and right panel of Figure [17] Thus, the distribution seen in 
Figure [17] (left panel) is due to an oversampling of the bright 
, young WD components that are more easily observed and 
can be found at farther distances. We will discuss the possible 
implications of this selection effect in more detail in Section 5 

4.5. Very close binary pairs 

In our sample of 1756 pairs, 1347 had individual expo- 
sures available in SDSS DR7 that were accessible to us from 
previous analyses (Wll). We found 37 objects that showed 
strong evidence of large changes in RV in small time inter- 
vals. The RVs and Lna/Lboi were calculated in the same 
way as was described in Section 3.1 and Section 3.3.5, re- 
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spectively. The objects' primary identifiers, median change 
in RV, median error in the change in radial velocity (con- 
strained using the method described in Section 3.1), and the 
median change in time in between exposures are reported in 
Table[4] As well as being interesting for studying the variabil- 
ity of WD+dM activity, these objects are excellent candidates 
for studying very cl ose binary pairs, mass ac cretion (e.g.. 



Hartig an et al.l |1995|). eclipsing binaries (e.g., iPvrzas et al.l 



rtiga 

20091 l2012t_ Parsons et al.l 120121). and pre-cataclysmic vari - 



if 

ables (e.g.. ISchreiber & Gansickell2003t iTappert et al1l2009h . 
We present these objects to the community as promising can- 
didates for follow-up observations. 

Two of the most interesting of the 37 objects can be found in 
Figure [18] For each pair, we plot the dM (triangles) and WD 
(squares) relative radial velocities vs. time elapsed since the 
first observation (hours). As discussed above, the WD RV cal- 
culations have larger uncertainties associated with them due 
to the broad WD absorption lines. Figure 1X81 also shows the 
Lna/Lboi vs. time elapsed in the bottom panel. Portions of 
each spectra containing the 8 183 A and 8194ANa I absorption 
doublet and the Ha line are included to provide further evi- 
dence for our radial velocity calculations. The vertical bars in- 
dicate the predicted location of the sodium doublet after being 
corrected to the measured system radial velocity. The fluxes 
of the spectra are normalized and offset for visualization pur- 
poses. 

The top right panel of Figure [18] shows the dM RV for 
SDSS J015225. 39-005808.6 changing by -200 km s" 1 over 
the course of an hour along with a correponding change 
in the WD RV by —100 km s _1 . Previous studies have 
identified SDSS J015225.39-005808.6 as a very close bi- 
nary system with a spectroscopically determin ed orbital pe- 
riod of 2.2h dNebot Gomez-Moran et al.l 1201 lT) . In the bot- 
tom plot of Figure [18] we see that the measured dM RV for 
SDSS J075919.42+321948.5 increases by -200 km s" 1 over 
an hour while the WD RV decreases by the same amount. 
We are likely seeing a significant component of the orbital 
motion in this pair in a short amount of time. Previous 
studies have also identified SDSS J075919.42+321948.5 as 
a very close binary system and candidate eclipsing binary 
dNebot Gomez-Moran et al.ll20Tll) . We suggest that the ob- 
jects found in Table [4] are promising cand idates for fol- 
lowup studies pertaining mass acc retion (e.g..lHartigan et ail 
19951) or eclisping binaries (e.g.. [Pvrzas et al.l 12 009. I20T2T 
Nebot Gomez-Moran et al. 2011tlParsons et al.ll2012l). 
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5. DISCUSSION 



We investigated how the magnetic activity of dMs in 
close WD+dM systems compared to that in the field pop- 
ulation (Wll). In Figure [10] we found the activity frac- 
tion in WD+dM pairs was higher than in isolated field 
dMs across all spectral types up to M7. We propose that 
the increased activity seen in WD+dM pairs is due to in- 
creased rotation i n the paired dMs. Many studies (e.g. , 
SkumanicrJ [19721: iPallavicini et alj|198il: iNoves et all 119841; 



Delfosse et al. f l 19981: iPizzolato et al J 120031; iBrowning et alj 
20101) have correlated stellar activity with rotation in main- 
sequence type stars to mid-to-late dM spectral types; the 
larger the rotation velocities, the more activity one expects. 
Thus, through disk-disruption, tidal effects, and angular mo- 
mentum exchange, we suspect that the rotation of the dM in 
a close WD+ dM binary sys t em is faster than that in a typ- 
ical field dM. iMeibom etal] ( 120071) argued that disk disrup- 
tion during the early stages of stellar dynamical evolution, 




Figure 17. WD cooling age vs. absolute Galactic height (left panel) and WD cooling age vs. distance (right panel). The methods used to calculate WD cooling 
ages (Gyr), absolute Galactic Heights (pc), and distances (pc) can be found in Section 3.2, 3.3.4, 3.3.3, respectively. The cluster of young WDs (WD cooling age 
~ 0.0) with a large spread in both Galactic height (left panel) and distance (right panel) indicates that our sample is biased towards WD+dM systems with young, 
bright WDs that we are capable of seeing at large distances. Galactic height and distance are roughly equivalent properties due to the SDSS footprint, which 
focuses on the Northern Galactic cap. We note that WD+dM pairs with distances < 300 pc span the entire range of WD cooling ages and are likely unaffected 
by the selection bias. 



Table 4 

Very close WD+dM candidate pairs 



SDSS ID 


Plate 


MJD 


Fiber 


#of 


ARV 


CTARV 


At 


s 












Exposures 


(km s" 1 ) 1 


(kms" 1 ) 2 


(hours) 3 


(AU) 4 


SDSS J001726.64- 


002451.2 


687 


52518 


153 


3 


42.1 


11.8 


0.375 


0.1564 


SDSS J002157.91- 


410331.6 


1913 


53321 


257 


8 


49.1 


14.7 


0.375 


0.0227 


SDSS J005245.12- 


-005337.2 


394 


51812 


96 


5 


42.1 


25.5 


0.281 


0.0115 


SDSS J015225.39- 


-005808.6 5 


1504 


52940 


83 


3 


219.0 


7.9 


0.375 


0.0758 


SDSS J023938.044 


-273654.1 


2444 


54082 


149 


10 


106.4 


38.4 


0.375 


0.0109 


SDSS J030138.244 


-050219.0 6 


2307 


53710 


140 


3 


64.9 


27.7 


0.469 


0.1476 


SDSS J030308.36-4 


-005444.1 7 


709 


52205 


524 


4 


330.9 


13.4 


0.469 


0.0088 


SDSS J073003.884 


-405450.1 


2683 


54153 


224 


6 


26.3 


7.5 


0.375 


0.0076 


SDSS J073455.924 


-410537.5 


2683 


54153 


507 


6 


31.1 


10.3 


0.375 


0.0282 


SDSS J075919.424 


-321948.5 


890 


52583 


631 


4 


51.1 


11.6 


0.375 


0.0955 



We calculated the change in radial velocities from one exposure to another (most observations were taken in succession) . ARV is the 
median of the changes. 

2 We report in this column the calculated the uncertainty in the reported ARV using each exposures individual radial velocity measurement 
and uncertainty. In order for the pair to be a very close WD+dM candidate ARV > <tarv- 

3 The time (hours) in between exposures for the reported ARV. 

4 Projected linear separation as defined in Section 3.3.4. Recall that the separations presented in this paper are upper limits. 

^ This system was identified as a close binary system in Nebot Gomez-Moran et al. 1 201 1 ) and has a measured spectroscopic orbital period 
of2.2h. 

^ This system was identified as a close binary system in Nebot Gomez-Moran et al. 1 20 1 1 ) and has a measured spectroscopic orbital period 
of 12.9h. 

7 This system was identified as a eclipsing binary in lPvrzas et all 120091) 



which prevents or limits angular momentum dissipation and 
spin-down of the dM, is more likely the culprit than is tidal 
synchronization. They found no correlation between stel- 
lar rotation and orbital period or eccentricities among stars 
in the young (~ 150 Myr) M35 open cluster. In addition, 
iMeibom et alj (|2P07!i) offered the suggestion that perhaps the 
observed faster rotation in close binaries could be due to the 
varying amounts of intrinsic angular momentum present dur- 
ing the formation and very early evolution of the closest bi- 
nary stars versus single stars. 

Previous studies (e.g., S06) also found that early-type dMs 
(< M4) with a close WD companion have higher activity frac- 
tions than isolated field dwarfs. Due to the small number of 
stars with spectral types > M5, S06 were unable to investi- 
gate the activity fractions of late-type WD+dMs in compar- 
ison with field dwarfs from W08. Our results for the early- 
type dMs agree with the results of S06 and our statistics are 
better with over 70,000 field dMs from Wl 1 as well as 1237 



WD+dMs with reliable activity properties from our sample 
(compared to 661 WD+dMs with activity properties from 
S06). The S06 study also lacked sufficient numbers to ex- 
plore the activity fractions in spectral types > M5 with any 
statistical significance. We have extended the activity fraction 
analysis in close WD+dM pairs out to spectral types ~ M7. 
We show with statistical significance that the trend of higher 
activity fractions for dMs in WD+dM pairs in comparison to 
field dMs extends at least to spectral types of M7 (Figure ITOb. 

An interesting feature worth noting from Figure |l0]is the 
convergence between the active fractions of binary and single 
dM populations towards later spectral types. Figure [TOlshows 
that the activity fraction in later spectral type dMs with close 
binary companions is not increasing as quickly as field dMs 
seen by Wl 1 . If the two populations truly converge at spectral 
types > M7, as Figure [To] suggests, then activity of late-type 
dMs in close WD+dM pairs increases by a smaller amount 
in comparison with field dMs. Due to the long lifetimes of 
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Figure 18. Two examples of the very close WD+dM pairs. Upper left pan- 
els: dM relative radial velocity vs. time elapsed since first observation. The 
black triangles are the dM radial velocity measurements while the squares are 
the WD radial velocity measurements. Lower left panel: L.H a /Lbol for each 
exposure. Middle panel: The 8183A and 8194A Na absorption lines from the 
dM for each exposure. The vertical bars indicate the radial velocity corrected 
location of the respective spectrum's Na I doublet lines, while the dashed 
lines indicate the rest wavelength of the Na absorption doublet (8 1 83 A and 
8194A ). Right panel: The dM spectrum centered on the Ha line. The fluxes 
in the middle and right panel are normalized and offset. 

M dwarfs and the fact that their activity state is correlated 
with their location in the Galaxy, it is more appropriate to an- 
alyze their activity fra c tion within the prop er Galactic context 
dKowalski et alJ 120091: Iffilton et alJ[2O10E W08. Wl 1). Fig- 
ure QT| confirms that dMs in close WD+dM pairs have an in- 
creased likelihood of being active. We assume that the main 
distinguishing property between dMs with close WD compan- 
ions and field dMs is that the former have more rapid rota- 
tion (magnetic locking, disk disruption, tidal effects, or an- 
gular momentum exchange) than we might expect. This im- 
plies that the mechanism driving magnetic field generation in 
late-type dMs may be less sensitive to rotation than the a-fi- 
dynamo thought to be responsible for magnetic fields in stars 
with masses greater than ~ 0.33 M©. Previous theoretical 
studies have developed approaches for generating magnetic 
fields in fully convective stars that are corr elated with rotation 
dChabrier & Kukerf2006MBrownind267j8l) . Observationally, a 
few studies have suggested that the correlation between rota- 



tion and activity in fully convective stars may not be as stron g 
as th eory predicts (I West & Basrill2009l IReiners et all 120 121) . 
Both IWest & Basril d2009l) and IReiners et al.l ( 120121) provide 
observations moderately rotating (> 3.5 km s _1 ), late-type 
dMs without any observable magnetic activity. While some of 
these previous results may be spurious due to low S/N spectra 
or other systematics, they offer a tantalizing possibility about 
the limited role of rotation in generating magnetic activity in 
late-type dMs. 

Figure [TT] also examines the time evolution of magnetic ac- 
tivity, where Galactic height was used as a proxy for the ap- 
proximate age of the system. We assume all of the stars are 
formed close to the mid-plane of the Galaxy and over their 
lifetimes they are systematically able to move further from 
the plane due to dynamical interactions. Figure [TT] shows 
that the activity fraction is significantly larger across all the 
Galactic height bins, indicating that not only are the early- 
type WD+dMs more likely to be active than the field popula- 
tion, but also have longer activity lifetimes (see W08). While 
the activity lifetimes in early-type WD+dMs may indeed be 
longer lived, Figure [16] (WD cooling ages as a function of 
spectral type for the active and inactive dMs) shows that there 
is a dearth of older active early-type dMs. The cooling age 
analysis in addition to the fact that not all early-type dM are 
active indicates that early-type dMs in close WD+dM bi- 
nary pairs still have finite active lifetimes (also seen in S06). 
This result is expected and agrees with the idea that dMs in 
WD+dM pairs are sped-up (with respect to field dMs) and 
still contain a solar-like a-il dynamo. However, for the late- 
type WD+dMs (> M5), the difference in the activity frac- 
tions becomes less obvious. In the M5 panel of Figure nil the 
WD+dM pairs trace the activity fraction of the field dMs very 
well, implying that the presence of a close WD companion no 
longer has an appreciable effect on the activity of dMs with 
spectral types > M5. Not only is the activity fraction consis- 
tent between the two populations at spectral type M5, but this 
similarity implies that the activity lifetimes are comparable 
between the two populations. This is corroborated by the WD 
cooling age analysis (Figure[T"6ll. where the active and inactive 
late-type WD+dM pairs appear to be drawn from the same 
parent distribution with higher confidence (82.3, 46.9, and 
100.0% for spectral type bins M4, M5, and M6+7) compared 
to the early-type WD+dM pairs (2.5, 1.7, and 2.5% for spec- 
tral type bines M0+M1, M2, and M3). With our assumption 
that dMs in close WD+dM pairs are rotating faster than the 
field population, this result suggests that rotation is not as im- 
portant in the activity state of late-type paired dMs. Evidence 
for this trend continues into the M6+M7 spectral types. Fig- 
ureQ~TJsuggests that whatever the driving mechanism for mag- 
netic field generation in late-type dMs, it may not be affected 
by a close binary companion, or even more dramatically, rota- 
tion. The absence of enhanced activity in late-type WD+dM 
systems can also be explained if the rotation-activity relation 
saturates at lower levels than in the early-type dMs. In this 
case, late-type dMs in close WD+dM pairs are rotating faster 
than their single counterparts. However, once the threshold 
is reached, increased r otation will not cause an increase in 
magn etic activity (e.g.. lPizzolato et aj]l2003t iBrowning et al.l 
l2010h . 

Figure [TT] shows that our sample may be biased towards 
younger systems, where both components are intrinsically 
brighter and easier to observe at greater distances. We there- 
fore, may be oversampling young WD+dM pairs, which (at 
least for early-type dMs) tend to be more active as seen in 
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Figures QT| and [16] Bins of higher Galactic Height (given 
the SDSS lines-of-sight are concentrated near the Northern 
Galactic Cap) will then be preferentially oversampled by ac- 
tive WD+dM stars simply because they are brighter and eas- 
ier to observe, tending to falsely inflate the activity fractions 
in the more distant Galactic Height bins. We examined this 
potential selection effect by remaking Figure QT|but on ly in- 
cluding WD+dM pairs with distances <300 pc, where we 
would be sensitive to WD+dM pairs of all WD cooling ages 
(Figure [19). We chose to cut on distances <300 pc because 
of the smooth sample distribution seen below this cutoff dis- 
tance (see Figure \Tt\ right panel). The distance cut reduced 
our sample to 892 WD+dM pairs, which has the consequence 
of enlarging the uncertainties in many bins. Figure [T9"1shows 
similar behavior to Figure [TT] namely increased activity frac- 
tion in the early-type dMs in WD+dM pairs in comparison to 
the field population (with the exception of one bin in M0+M1 
which is likely due to small number statistics) and similar ac- 
tivity fractions between the late-type WD+dM and field dM 
populations. The trend of decreasing activity with Galac- 
tic height is somewhat lost in the M2 and M3 spectral type 
bins due to the small numbers and the fact we are sampling 
a smaller range of Galactic heights. We conclude that hav- 
ing a preferential selection bias towards brighter and therefore 
younger systems at large distances does not make a qualitative 
difference in our results. 
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Figure 19. Fraction of active stars as a function of absolute Galactic height 
(pc) for each spectral type for stars with distances < 300 pc. The distance 
selection was introduced to test the severity of the selection bias identified 
by Figure 1171 In general, the lack of distant, young systems does not affect 
the activity fractions and we conclude that the selection bias does not signifi- 
cantly change our results. 

Figures[lO]and[TT]explored the behavior of the activity state 
(and activity lifetimes) of dMs with close WD companions. 
We found evidence that the presence of close companions in- 
creases the likelihood that a star is active and extends the life- 
times of early-type dMs. However, for late-type dMs, close 
companions did not have much of an effect on the activity 
and activity lifetimes of late-type dMs. We used Figure [T2"1to 
make an analogous study, except exploring activity strength 
(via Lna/Lboi) of close WD+dM pairs in comparison to 
the field dM population. In Figure [12] the median activity 
strength of the active WD+dM population is slightly higher 
across all spectral types in comparison to the field dM pop- 



ulation. Also, Lna/Lboi is shown to remain relatively con- 
stant at approximately Lna/Lboi = 10 -3- ' from spectral types 
M0- M4; a trend that is consistent with the f ield dM popula- 
tion (iBurgasser et al.ll2002t iGizis et al.ll2002l W04). Towards 
the later spectral types (> M5), there is a decreas e in activity 
in both populations (also seen in previous studies. IGizis et all 
l200l IWest etani200l 120061 Wll). In the M6+M7 spectral 
type bin we see a distinctly larger difference in Lna/Lboi be- 
tween the two populations. The question is whether or not the 
M6+M7 result is real or simply due to a small number of stars 
in the M6+M7 spectral type bin. If it is a real result, then it 
could possibly provide some interesting insight into the mag- 
netic field generation in stars with fully convective envelopes. 
Figures [tOland [TTIshow that the activity state in later spectral 
types does not seem to be affected by having a close binary 
companion, however, the strength of activity might be ampli- 
fied by the WD. If having a close companion imparts a larger 
rotation to the dM, then we could be seeing that rotation in 
later spectral types influences the strength in activity, but has 
less of an effect on the likelihood of activity. In other words, 
rotation may play a limited role establishing magnetic activity 
but may have an effect on the strength of the chromospheric 
heating. 

Figure Qj] puts strength for each spectral type bin in 
the proper Galactic context. Figure [13] shows that the 
LHa/Lboi for M0+M1 and M2 dwarfs appears to be mostly 
flat for both the WD+dM and field populations for most 
Galactic heights, but that the field dMs at spectral types of 
M3 and M4 show a slight decrease in activity with increas- 
ing Galactic height. There is no clear trend of decreasing 
LHa/Lboi with Galactic height seen in spectral types M3-M4 
for the WD+dM population. At spectral types of M5 there 
is a large decrease in overall activity strength in both pop- 
ulations, and a distinct decline in LHa/Lboi with increasing 
Galactic height in the field population, a trend that is indica- 
tive of activity strength decreasing with age. There is no ap- 
parent trend in the WD+dM population that shows a decrease 
in LHa/Lboi with Galactic height as seen in the field dMs. 
While there are only a few objects in the M6+M7 bin, it is cu- 
rious that there is a fairly large difference in activity strength 
between the two populations, as was identified in Figure Q~2] 
We reiterate that this could be evidence for close companions 
affecting the activity strength in dMs but not increasing the 
likelihood of activity or the lifetime of activity (Figure LLD. 

We have shown convincing evidence that close companions 
have an effect on both the activity state and activity strength in 
dMs. As was discussed in Section 3.3.2, we have an approxi- 
mate measure for the separation for which we can investigate 
the effect of separation on the activity of dMs. Figure [l4lis 
presented in the same way as Figure [10] where the activity 
fraction for the WD+dM population is shown as a function of 
spectral type. The active WD+dM population is further split 
into three separate separation bins, those less than 0.1 AU, 
0.1-1 AU, and 1-100 AU. These separations bins were chosen 
to ensure good statistics in each spectral type bin. In com- 
paring Figure [l4]to Figure [10] the activity fraction across the 
spectral types follows the same trend for each separation bin 
as it did for the total population: low activity fractions in early 
spectral types and increasing activity fractions toward the later 
spectral types. We find that the closer separated pairs are more 
active than the wider separation pairs, while the widest pairs 
are still more active (and likely to be active) than their field 
counterparts. We also report an interesting trend that is seen 
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at the later spectral types, where the active fractions for three 
different separation populations, as well as the field popula- 
tion appear to converge. Together with other results in this 
paper, we suggest that the mechanism driving magnetic fields 
in late-type dMs becomes decreasingly sensitive to the pres- 
ence of a close companion, which we have interpreted as an 
insensitivity to rotation. However, Figure [l5]shows that while 
every spectral type bin shows signs of decreasing activity with 
increased separation, a linear least-squares fit is only at most 
50% significant over a flat horizontal line drawn through the 
data for each spectral type. Figure [14] shows that the activity 
state is highly dependent on the separation of the WD+dM bi- 
nary system while Figure[l5]shows that evidence for a correla- 
tion between activity strength and separation of the WD+dM 
binary system is inconclusive. As a result of our study and 
compilation of the WD+dM sample we have found a number 
of objects that have shown convincing signs of being in a very 
tight binary system. We have identified 37 objects in Table |4] 
that we feel may be valuable to the astronomical community 
in studying dynamics of close binary systems, mass accretion, 
eclipsing binaries, and pre-cataclysmic variables. 

6. SUMMARY 

Using close WD+dM pairs, we have investigated the mag- 
netic activity properties of M dwarfs in WD+dM pairs. We 
used SDSS DR8 spectra of unresolved binaries and developed 
an iterative procedure to separate the two spectral components 
in 1756 high S/N WD+dM pairs. We used the Ha line as 
an indicator of magnetic activity and both absolute Galactic 
height (pc) and WD cooling models to examine the activity 
properties as a function of age. Both the activity state and the 
activity strength (via Lna/Lboi ) as a function of separation 
were investigated using approximate binary separations found 
from the dM and WD component radial velocities. Below we 
summarize the findings of this paper. 

1. We built upon previou s clos e WD+dM investigations 
(iRavmond et alJ 12003b S06; iRebassa-Mansergas et ail 

and compiled a sample of 1756 high S/N close 
WD+dM pairs. 

2. We calculated numerous parameters for our sample in- 
cluding: dM spectral type, WD effective temperature, 
WD surface gravity, WD and dM relative radial veloci- 
ties (km s _1 ), the system radial velocity through space 
(km s _1 ), galactic height above the plane (pc), approx- 
imate distance to the system (pc), approximate binary 
separation (AU), period of binary (days), activity (as 
traced by Ha), and Lna/Lboi, and WD ages from cool- 
ing models. 

3. We confirmed that paired dMs are more active than field 
dwarfs for early dM spectral types but then become 
comparable in their activity fractions at later spectral 
types. This gives further evidence to the notion that 
the mechanism that generates magnetic fields in late- 
type dM spectra may be different than the mechanism 
in early-type stars, and that it may be relatively unaf- 
fected by a close companion. 

4. The activity lifetimes of early-type paired dMs is longer 
than field dMs, but still shows sign of decreasing with 
increasing height, implying a finite lifetime. However, 
late-type paired dMs have lifetimes consistent with the 
lifetimes found in isolated dMs. This is reiterated when 



using WD cooling ages and finding that early-type ac- 
tive dMs in WD+dM pairs are on average younger and 
part of a distinct population in comparison to inactive 
early-type dMs in WD+dM pairs. In contrast, the late- 
type active dMs in WD+dM pairs span a broader range 
in age and are indistinct from the late-type inactive dMs 
in WD+dM pairs. 

5. We found that the median activity strength in WD+dM 
pairs is larger across every spectral type when compar- 
ing to the field dM population. When analyzing ac- 
tivity strength in the proper Galactic context we find 
that spectral types M0-M2 Lna/Lboi remains con- 
stant for both populations; M3-M4 field dMs exhibit 
a decreasing Lna/Lboi with increasing galactic height 
while WD+dMs do not exhibit a decreasing trend and 
are still elevated in Lna/Lboi in comparison with the 
field population; M5 field dMs again exhibit a de- 
creasing Lna/Lboi with Galactic height, again not seen 
in WD+dMs, M6+M7 field dMs show a significantly 
lower Lna/Lboi than in earlier spectral types while 
M6+M7 did not decrease in Lna/Lboi as significantly 
as the field population. 

6. The closer the separation between WD+dM (and pre- 
sumably any X+dM system) the more likely the dM is 
to be active. 

7. The activity strength in WD+dMs shows no conclusive 
evidence for being dependent upon the proximity of the 
binary pair. 

8. We present 37 candidate very close pairs exhibiting sig- 
nificant radial velocity changes on hour timescales to 
which we hope will be useful for the community for 
follow up studies. 

In investigating the activity properties of M dwarfs with 
close WD companions we have found some interesting trends 
occurring near the transition at which dMs become fully 
convective. Spectral types ranging from M0-M7 have been 
shown to exhibit higher activity fractions than the field dM 
population (Figure [TOl. However, a possible convergence be- 
tween the WD+dM and field dM population is seen at later 
spectral types. When put in the proper Galactic context, we 
show that early-type WD+dMs still exhibit higher activity 
fractions as well as extended (albeit finite) activity lifetimes 
when compared to the field population. Beginning at spec- 
tral type of M5, we find that the activity behavior in WD+dM 
becomes very similar to that of the field population; both in 
activity fractions and lifetimes. This change conveniently oc- 
curs around the proposed ^M4 transition to fully convective 
envelopes. We argue that the primary difference between the 
WD+dM and field population is that the rotation of the dM 
is being sped up by the presence of a close WD compan- 
ion; most likely through the destruction disk and thus reduc- 
ing magnetic brakin g in the dM and strengthening rotation 
iMeibom et al.l (120071) . The decrease in activity with the sup- 
posed increased rotation in late-type WD+dM implies that 
the magnetic field generation mechanism is not as sensitive to 
rotation as in early-type M dwarfs. Almost in contradiction, 
we find that the activity strength does seem to increase in the 
presence of a WD companion in late-type WD+dM. We have 
also shown that Lna/Lboi increases with the proximity of the 
binary companion. With our results we suggest that rotation 
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in late-type dMs has a bimodal effect; rotation may play less 
of a role in establishing magnetic activity but does have an 
effect on the strength of the magnetic activity. 

There are, however, alternative solutions. The first one be- 
ing that rotation is important for activity and the proximity of 
a close binary companion is not a reliable indicator of the dM 
being "sped" up in comparison to a field dM. Another solution 
is that at later spectral types there is a lower rotation thresh- 
old that sets the activity and most of the late-type stars in both 
the populations are above th at threshold. This so lution is put 
into question by the study by lWest & Basril (|2009), which pre- 
sented field M7 dMs with significant rotation but no measured 
activity. Also, Ha emission can be induced on the dM by the 
irrad iation of a hot WD companion in close binary systems 
(e.g. lTappertetai1l20lU iParsons et ail 120 121) . However, the 
effects of the irradiation o nly becomes signific ant for systems 
with WD T eff > 45, 000K (iTappert et alJ201 lh . Only 78 (4%) 
of our sample have T e s > 45,000 K with the majority having 
T eff between 10000K and 25000K. Hence, Ha emission in 
dM induced by WD irradiation is expected to be negligible 
for 96% of our sample; we expect it to have little effect on the 
results presented in this paper. 

In the future, we hope to increase the sample size of 
WD+dM pairs with spectral types > M6 and extending this 
analysis to lower masses in addition to improving the over- 
all statistics. Despite the current dearth of late-type dMs in 
WD+dMs companions, we have added to previous studies 
that find likely differences in the dynamo mechanism across 
the fully convective boundary. We believe that close WD+dM 
pairs provide some unique insight into the behavior of mag- 
netic fields across the fully-convective stellar mass regime and 
will aid in constraining the models for generation and main- 
taining magnetic fields in fully convective M dwarfs. 
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